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INTRODUZIONE 


L'Italia, negli ultimi anni, è stata colpita da intense attività sismiche che hanno coinvolto 
soprattutto il centro della penisola, producendo molte vittime e causando enormi danni 
al patrimonio artistico e architettonico. A essere oggetto di crolli sono stati soprattutto 
edifici in muratura caratterizzati da forti irregolarità geometriche e realizzati con materiali 
a volte scadenti, che mal si comportano sotto l’effetto delle azioni sismiche. Da questo 
punto di vista, gli edifici maggiormente colpiti sono le chiese, che per loro natura non 
riescono a sviluppare un comportamento scatolare generando quadri fessurativi 


importanti nonché crolli nei casi peggiori. 


Il problema della sicurezza e della preservazione di questi importanti edifici storici rimane 
una questione aperta e, pertanto, la valutazione della loro risposta ad azioni dinamiche 


risulta di essenziale importanza, specialmente nelle aree con rischio sismico elevato. 


Negli ultimi anni, attraverso lo sviluppo di tecnologie software e hardware sempre più 
avanzate nell’ambito dell’Ingegneria Civile, sono stati sviluppati sistemi di Structural 
Health Monitoring (SHM), cioè di valutazione dello stato di salute strutturale in tempo 
reale, con l’obiettivo di cogliere rapidamente eventuali variazioni nel comportamento e 
nelle principali caratteristiche della struttura al fine di valutarne costantemente lo stato 
di salute e anche di fungere da utile supporto nell’individuazione del danno. Questo 
monitoraggio nel tempo può basarsi anche sul tracciamento del comportamento 
dinamico della struttura nel tempo, quindi sia in condizioni normali di esercizio (ossia 
struttura sottoposta a vibrazioni di tipo ambientale) che a seguito di eventi eccezionali 


(come nel caso dei terremoti). 


Il presente lavoro di tesi ha come obiettivo l’analisi e lo studio dei risultati preliminari di 
un sistema di monitoraggio dinamico della salute strutturale messo a punto nella chiesa 
di Santa Maria in Via in Camerino (MC), gravemente danneggiata in seguito al sisma del 


Centro Italia del 2016. Inoltre, illavoro ha visto anche l’implementazione di una procedura 
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automatica per l'individuazione dei parametri monitorati sulla base delle registrazioni 


acquisite, parametri sui quali si basa il monitoraggio. 


Da un punto di vista operativo, nel primo capitolo, dopo una breve introduzione sui 
problemi del patrimonio architettonico italiano in termini di vulnerabilità sismica, si 
riporta una panoramica di diversi casi studio di chiese condotti soprattutto in Italia, ma 
anche in Europa. In primo luogo, sono stati esaminati articoli di ricerca condotti in edifici 
lesionati dal terremoto dell’Aquila del 2009 come la Basilica di Santa Maria in Collemaggio, 
e le Chiese di San Giuliano, di Santa Maria del Suffragio, di San Pietro a Coppito e la Basilica 
di San Nicola da Tolentino danneggiata dal Sisma del 2016. Nella seconda parte sono state 
prese in considerazione strutture di notevole rilevanza architettonica e culturale, come il 
Duomo di Milano, il Santuario di Vicoforte a Cuneo, il Battistero della chiesa di San 
Giovanni, la Cattedrale di Colonia, la Chiesa di San Torcato a Guimaraes, la Chiesa del 
Monastero dei Girolamini a Lisbona e la Cattedrale DI Santa Maria Assunta a Reggio Emilia 
focalizzando l’attenzione sul tipo di indagini effettuate e sui risultati ottenuti. Vale la pena 
ricordare che essendo questa metodologia in via di sviluppo, non sono molti gli studi 
presenti in letteratura e la ricerca condotta nella chiesa di Santa Maria in Via può costituire 


la base per successivi sviluppi della tecnologia SHM. 


Nel secondo capitolo si è illustrato il caso studio della chiesa di Santa Maria in Via 
ripercorrendone le fasi costruttive con particolare attenzione agli eventi sismici che hanno 
interessato l’edificio nel corso della sua esistenza. Quindi si è descritto il danneggiamento 
conseguente allo sciame sismico del Centro Italia del 2016 che oltre a un ingente quadro 
fessurativo ha prodotto anche il crollo del campanile nel retro della chiesa e ha innescato 
un meccanismo di ribaltamento al piede della facciata. Inoltre, è stato descritto con 
attenzione il presidio di sicurezza in acciaio messo in opera per contrastare il cinematismo 


attivato. 


Nel terzo capitolo è stato descritto nel dettaglio il sistema di Structural Health Monitoring 


implementato per la chiesa in questione. Inizialmente è stata delienata l’identificazione 
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dinamica globale eseguita attraverso l’analisi modale operativa (OMA) che ha portato 
principalmente all’individuazione di sei forme modali. Si è così predisposto un sistema di 
monitoraggio strutturale al fine di monitorare i modi di vibrare individuati 
preliminarmente. Sono state descritte scrupolosamente sia la strumentazione adottata 
che le operazioni svolte sulle registrazioni acquisite. Si è illustrato, altresì, l'algoritmo SSI- 
PC adottato nell’identificazione dinamica e particolare attenzione è stata dedicata 
all’identificazione automatica dei parametri modali basata sulla tecnica del clustering 


gerarchico agglomerativo che svolge un ruolo cruciale nei sistemi di SHM. 


In conclusione, nell’ultimo capitolo vengono riportati i risultati del monitoraggio eseguito 
dal 16 novembre 2020 al 6 gennaio 2021 e quanto emerso dall’identificazione automatica 
dei segnali acquisiti. Si riporta infine un’analisi statistica descrittiva delle forme modali 
identificate e sono state analizzate le frequenze in relazione ai parametri ambientali 


acquisiti. 
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1.LO STATO DELL'ARTE E LA RICERCA DI 
LETTERATURA 


Il patrimonio architettonico italiano rappresenta una fonte di inestimabile ricchezza e 
bellezza che deve essere protetta. Negli ultimi anni, importanti eventi sismici hanno 
colpito la penisola, tra i più importanti si ricordano: Umbria Marche del 1997, San Giuliano 
di Puglia nel 2002, L'Aquila 2009, Finale Emila del 2012 e il terremoto del Centro Italia del 
2016. Alla luce di questo scenario che ha generato molte vittime e distrutto molte realtà 
locali, è emerso il problema della sicurezza degli edifici esistenti sia in muratura che in 
calcestruzzo armato. Questo è dovuto in particolar modo all'elevato rischio sismico del 
territorio italiano che viene quantificato in termini di pericolosità, esposizione e 
vulnerabilità. Proprio su quest’ultimo aspetto si deve agire al fine di accrescere il livello di 


sicurezza dell’edificato italiano. 


Per questi motivi, è stata condotta un'intensa attività di ricerca che ha riguardato 
principalmente strutture storiche come le chiese, che sono molto vulnerabili da un punto 
di vista sismico e che d’altra parte rappresentano un’inestimabile ricchezza che non può 
svanire nel tempo. Da questo punto di vista il monitoraggio dinamico dello stato di salute, 
SHM, può offrire una soluzione ottimale per il controllo in fase di esercizio, mentre su 
quegli edifici danneggiati dal sisma può essere uno strumento utile attraverso il quale 
pianificare eventuali interventi di riabilitazione strutturale. Sono quindi riportati gli aspetti 


principali di alcuni articoli di ricerca condotti in Italia e in Europa. 


Li: Il rischio sismico 


Per rischio sismico si intende la misura dei danni attesi in un dato intervallo di tempo. Al 
fine di quantificare ciò si deve considerare la combinazione di tre elementi: la pericolosità 


sismica, la vulnerabilità e l'esposizione. 


La pericolosità sismica rappresenta la probabilità che si verifichi nel territorio un evento 
sismico di una certa magnitudo in un certo intervallo di tempo. Si distingue in pericolosità 
sismica di base, che identifica il moto sismico di riferimento per il sito considerato, nelle 
condizioni ideali di terreno rigido e superficie topografica orizzontale, e pericolosità 


sismica locale che identifica il moto sismico di riferimento per il sito considerato. 


La vulnerabilità sismica è una misura della propensione di una struttura a subire un 


danno di un determinato livello, a fronte di un evento sismico di una data intensità. 


La terza componente del rischio è l'esposizione che viene intesa come la quantità e la 
qualità di elementi antropici esposti all'evento sismico, inclusa la popolazione (Corradi, 


Salvucci, & Vitale., 2015) 


In Italia questi fattori sono medio-alti in quanto è una zona fortemente sismica e sono 
presenti molti edifici in muratura che resistono poco alle azioni di trazione trasmesse dalle 


azioni sismiche. 


Si può ridurre il rischio sismico attraverso una riduzione della vulnerabilità in quanto non 
è possibile intervenire né sull’esposizione né sulla pericolosità. Ciò può essere raggiunto 
adeguando le costruzioni esistenti ai nuovi standard di sicurezza, ottimizzando la 
manutenzione, monitorando gli edifici esistenti e migliorando anche i servizi di emergenza 


locali nel caso dovessero verificarsi degli eventi. 


LZ; La ricerca di letteratura 


È stata condotta un’ampia ricerca bibliografica in quanto le prove dinamiche inizialmente 
furono applicate principalmente allo studio dei ponti e negli edifici di nuova costruzione, 
mentre negli ultimi anni grazie alla loro versatilità ed economicità sono state impiegate 
anche nello studio di strutture esistenti, edifici lesionati dal sisma e siti di importante 


valenza storico culturale. 


In alcuni di questi studi sono stati effettuati monitoraggi statici, in altri sono state 
realizzate delle prove di identificazione dinamica globale attraverso l’analisi modale 
operativa, mentre solo in alcuni di essi è stato applicato un sistema per il monitoraggio 


dello stato di salute SHM. 


In virtù di queste considerazioni, prima sono stati presi in esame i casi studio di edifici 


lesionati dal sisma come: 


e La Chiesa di Santa Maria in Collemaggio a L'Aquila; 


e La Chiesa di San Giuliano a Poggio Picenze; 


e La Chiesa di Santa Maria del Suffragio a L'Aquila 


e La Chiesa di San Pietro a Coppito a L'Aquila; 


e La Basilica di San Nicola da Tolentino; 


Nella seconda parte sono stati analizzati i casi studio riguardanti edifici di pregio 


architettonico e culturale o con specifiche problematiche costruttive: 


e Il Duomodi Milano; 


e |lSantuario di Vicoforte a Cuneo; 


e ll Battistero della Chiesa di San Giovanni a Firenze; 


e La Cattedrale di Colonia in Germania; 


e La Chiesa di San Torcato a Guimaràes in Portogallo; 


e La Chiesa del monastero dei Girolamini a Lisbona in Portogallo; 


e La Cattedrale di Santa Maria Assunta a Reggio Emilia. 


In Figura 1 si riporta una mappa dei casi studio presi in considerazione. 
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Figura 1: Mappa casi studio analizzati. 


LZ Prove dinamiche applicate a edifici lesionati dal sisma 


Nel presente paragrafo vengono analizzati nel dettaglio gli articoli di ricerca inerenti a 
prove dinamiche applicate a edifici lesionati dal sisma dell’Aquila del 2009 e al sisma del 
Centro Italia del 2016. Per ogni caso studio viene prima fatto un accenno alla sua storia, 


poi viene descritta la sperimentazione effettuata quindi i risultati e le conclusioni emerse. 


1.2:1.1. La Basilica di Santa Maria in Collemaggio a L’Aquila 


Figura 2: La Chiesa di Santa Maria in Collemaggio. 


Cenni storici 


La Basilica di Santa Maria di Collemaggio a l'Aquila, rappresenta un esempio di 
architettura romanica in Abruzzo. È stata costruita nel 1287 su una collina all'ingresso 
della città di L'Aquila. Nei secoli la struttura ha subito numerose modifiche e continue 
ristrutturazioni anche in seguito ai tanti terremoti che hanno colpito la zona. Nei primi 
anni 2000 sono stati eseguiti degli interventi al fine di ridurre la vulnerabilità sismica delle 
pareti longitudinali mediante iniezioni di malta cementizia ed è stato installato un sistema 
di controventatura in acciaio sotto la copertura della navata centrale al fine di limitare il 


trasferimento di azioni sulle pareti convergenti. 


La struttura è stata gravemente danneggiata in seguito al sisma del 6 aprile 2009. | danni 
principali sono stati osservati nella zona del transetto, interessata da un meccanismo di 
implosione di tutte le strutture che lo compongono: arco trionfale, cupola, le strutture di 
copertura e i pilastri. Nei primi mesi successivi al sisma, la basilica venne messa in 
sicurezza attraverso l'installazione di una copertura provvisoria. In seguito ai lavori di 


restauro venne riaperta al pubblico il 13 dicembre 2017. 


Descrizione della sperimentazione effettuata 


Dopo il terremoto è stato sviluppato e installato un sistema SHM con sensori posizionati 


sia internamente che esternamente. 
Gli obiettivi principali di questa ricerca sono stati: 
e  L’investigazione in merito alle possibili cause del collasso; 


e |lmonitoraggio della struttura metallica e di altri interventi di messa in sicurezza 


provvisori; 
e L'individuazione preventiva di un peggioramento eventuale dello stato di danno; 


e  L’esecuzione di un’analisi a lungo termine della risposta dinamica della struttura 
includendo gli effetti ambientali stagionali e le modifiche di tali parametri dopo il 


restauro. 


Tabella 1: Strumentazione adottata per il sistema di monitoraggio dinamico. 


16 ACCELEROMETRI TRIASSIALI (MEMS-LIS334ALH) 
8 ESTENSIMETRI 
3 INCLINOMETRI 
SCHEDE CONVERSIONE SEGNALE A/D 16 BIT (QuikFilter QF4A512) 
WIRELESS PLATFORM (MEMSIC IMOTE2 MOTE) 
SENSORE TEMPERATURA E UMIDITA' (SHT11) 
SENSORE DI LUMINOSITA' (TAOS 2561) 


La strumentazione adottata è stata posizionata sia all’interno che all’esterno della 
struttura, le caratteristiche sono ripotate in Tabella 1 mentre la disposizione planimetrica 


viene riportata in Figura 3. 


Layout of the complete structural 
health monitoring system 


Figura 3: Disposizione dei sensori: 16 accelerometri di cu 15 in elevazione rappresentati da un cerchio rosso e 
uno alla base individuato da un cerchio verde; 8 estensimetri denotati con un quadrato blu; 3 inclinometri a 
muro denotati con un triangolo arancione; gateway centrale rappresentato da una stella. 


Il sistema di monitoraggio è stato effettuato con sensori Wireless ed è stato in funzione 


per tre anni, segnatamente dal 2011 al 2014. 


Il lavoro è stato condotto in tre fasi principali: calibrazione in laboratorio, validazione 


iniziale in sito dei risultati, gestione e miglioramento della rete da remoto. 


La prima acquisizione accelerometrica nella basilica sotto sollecitazioni ambientale ha 
mostrato che la struttura in muratura, supportata dall’impalcatura metallica, era soggetta 
a bassissimi valori di ampiezze di vibrazione, di conseguenza era necessario un 
significativo input di energia per eccitare sufficientemente tutte le principali componenti 
modali. Pertanto, è stato condotto uno studio sul comportamento della struttura sotto 
l’effetto di otto eventi sismici che si sono verificati durante il periodo in cui il monitoraggio 


è stato in funzione. (Potenza, et al., 2015) 


| dati registrati sono stati processati mediante gli algoritmi SSI-COV e SSI-Data. 


Con i dati del sistema di monitoraggio è stato calibrato un modello agli elementi finiti della 


basilica che tenesse conto dello stato di danno attuale. 


Il processo di updating del modello è stato condotto variando sia le caratteristiche dei 
materiali che il grado di vincolo interno tra pareti ed il muro. In Figura 4 sono riportate le 


caratteristiche del modello effettuato. 
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Figura 4: Modello FEM della basilica. 


Risultati e conclusioni 
Dallo studio condotto, è stato possibile evidenziare i seguenti aspetti. 


In primo luogo, le misure dei dati registrati negli otto eventi sismici hanno dimostrato 
l'efficacia dell'algoritmo SSI nella valutazione e nell’identificazione dei parametri modali 
portando infatti all’identificazione di tre forme modali. Affinché ciò sia possibile è 


necessaria una scelta di sensori altamente performanti. 


Il trattamento attento dei dati acquisiti da diverse fonti di eccitazione ha permesso di 
eseguire l'aggiornamento di un modello agli elementi finiti della basilica nelle condizioni 
di danno attuali. Dall'analisi dei risultati riportati in Tabella 2 è possibile notare una 
discreta corrispondenza tra le misure delle frequenze identificate e quelle relative al 


modello numerico. (Potenza, et al., 2015) 


In Figura 5 si riportano le forme modali emerse dal modello agli elementi finiti. 


Tabella 2: Updating del modello: frequenze e MAC. 


Modes Frequencies MAC 

Identified Numerical D % Reduced Expanded 
1 0.9780 0.9500 2.95 0.9885  0.2081 0.5100  0.9929 0.0078  0.0010 
2 1.3970 1.5600 —10.45 0.4118 0.9727 0.0555 0.0079 0.9706 0.0209 
3 1.9100 2.0100 —4.98  0.1327  0.0940 0.7840 0.0343 0.0240 (0.9447 


Second mode 


Third mode 


(b) (ec) 


(d) 


Figura 5: Prime tre forme modali con viste in sezione e in pianta. 


1.2.1.2. La Chiesa di San Giuliano a Poggio Picenze (Aq) 


Figura 6: La Chiesa di San Giuliano a Poggio Picenze (Aq). 


Cenni storici 


La Chiesa di San Giuliano a Poggio Picenze in provincia di L'Aquila è stata costruita attorno 
al XV secolo sulle rovine di un ospedale rimasto operativo fino al 1447. Nella facciata è 
riportata la data del 1589 che probabilmente indica il termine della costruzione. Presenta 
sobria architettura, con il fronte nitido e con lesene in pietra locale, l'interno è a navata 


unica, con tetto a capriate lignee. 


La chiesa in seguito al terremoto dell'Aquila del 2009 ha subito notevoli danni ed è stata 
interessata da un pesante quadro fessurativo. Sulla facciata era presente un meccanismo 
di ribaltamento fuori dal piano che ha prodotto alcune fessure nelle zone di intersezione 
con le pareti longitudinali. Inoltre, sempre in questa zona erano presenti anche delle 
fessure a taglio e delle crepe verticali. È stata coinvolta anche la parte longitudinale della 
chiesa, in dettaglio è interessato il muro perimetrale sul lato destro che è soggetto a un 


ribaltamento locale causato dalla spinta di una struttura adiacente. 


Descrizione della sperimentazione effettuata 


Sono stati eseguiti dei test di vibrazione ambientale (AVT) per la definizione dello stato di 
salute della chiesa. La strumentazione utilizzata per la realizzazione di queste indagini 


viene riportata Tabella 3 e in Figura 7. 


Tabella 3: Caratteristiche principali della strumentazione adottata per effettuare i test di vibrazione 
ambientale (AVT). 


3 SISMOMETRI (SS-1 Ranger, Kinemetrics) 


1 CONDIZIONATORE DI SEGNALE A 4 CANALI 
1 SISTEMA DI ACQUISIZIONE DATI AD ALTA VELOCITA' 
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Kinemetrics product Four Channel Signal Conditioner High-speed data acquisition system 


Figura 7: Strumentazione adottata. 


La procedura dell'analisi di vibrazione ambientale ha previsto una fase iniziale di 
calibrazione durante la quale tutti i sensori sono stati posizionati nella medesima 
posizione e direzione. Il segnale quindi è stato registrato simultaneamente ed è stato 


ottenuto uno spettro di Fourier. 


Dopo questa fase i sensori sono stati posizionati a differenti altezze e in differenti punti di 
misura per ottenere le forme modali. Nel caso in esame i segnali registrati in diversi punti 
della struttura sono normalizzati rispetto ad un punto specifico denominato “reference 
point”, indicato in blu in Figura 8, dove è stata effettuata in precedenza la calibrazione. Il 
tempo di acquisizione della registrazione è stato di mille secondi. (Di Lorenzo, Formisano, 


Krstevska, & Landolfo, 2019) 


Figura 8: Disposizione dei sismometri nell’edificio. 
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I segnali acquisiti sono stati processati tramite il software ARTeMIS e processati mediante 
l'algoritmo FDD. In totale sono state effettuate misure in undici diversi punti della chiesa 


riportati in verde in Figura 8. 


I parametri sperimentali così ottenuti sono stati utilizzati per calibrare un modello agli 


elementi finiti della struttura mediante il software ABACUS come illustrato in Figura 9. 


Longitudinal mode with local effects: 
f=6.70 Hz 


Longitudinal mode: f =5.30 Hz 


Torsion mode: f =7.20 Hz 


Figura 9: Forme modali della struttura ottenute mediante il software Abacus. 


Risultati e conclusioni 


Le prove di vibrazione ambientale condotte nella Chiesa di San Giuliano hanno avuto 
come obiettivo quello di ottenere i parametri modali intesi come frequenze naturali, 


forme modali e rapporti di smorzamento, in seguito al terremoto dell’Aquila del 2009. 


In Figura 10, si riporta il Peak Picking delle frequenze dominanti. In corrispondenza del 
primo valore singolare, si notano tre picchi sostanzialmente che vengono associati a dei 
modi globali, mentre i picchi evidenziati nel secondo valore singolare possono essere 
associati a dei modi locali. In Tabella 4 si riportano le frequenze principali e i rapporti di 


smorzamento connessi alle forme modali individuate. 
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Figura 10: Peak Picking delle frequenze dominanti. 


Tabella 4: Frequenze dominanti e rapporti di smorzamento. 


Mode shape Dominant frequency (Hz) Damping 
coefficient 
%) 
Longitudinal 5.37 2.2 
Longitudinal 6.64 2.2 
Torsion 7.03 1,9 


Grazie alla campagna di indagine realizzata, si sono potute rilevare tre frequenze modali 
di cui le prime due corrispondenti a dei modi di vibrare flessionali mentre la terza è 
prettamente torsionale a causa anche del confinamento offerto dalle strutture adiacenti 
alla chiesa. Inoltre, vale la pena notare che, in virtù dell'importante quadro fessurativo gli 
spettri spesso contengono delle frequenze relative ad alcune parti di struttura 


danneggiate che di fatto complicano l’identificazione dinamica globale. 


Emerge infine, una buona corrispondenza tra i dati ottenuti sperimentalmente e i 


parametri numerici derivanti dal modello agli elementi finiti. 


Tabella 5: Confronto tra le frequenze sperimentali e numeriche. 


Mode Experimental frequency Numerical 
(Hz) frequency 
(Hz) 
Longitudinal 5.37 5.30 
Longitudinal 6.64 6.70 
Torsion 7.03 7.20 
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L:2:1°3; La Chiesa di Santa Maria del Suffragio a L'Aquila 
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Figura 11: Chiesa di Santa Maria del Suffragio a L'Aquila. 


Cenni storici 


La Chiesa di Santa Maria del Suffragio venne fatta edificare nel 1713 in suffragio delle 
vittime del terremoto del 1703 e terminata nel 1803 con la realizzazione della cupola 
attribuita a Giuseppe Valadier. Il terremoto del 6 aprile del 2009 ha generato un nuovo 
collasso della chiesa. Il quadro fessurativo era caratterizzato da fessure a taglio e a 
compressione generate dal movimento a pendolo della struttura e dal contemporaneo 
crollo di parte di essa. | pilastri più danneggiati erano quelli dalla parte del transetto con 
fessure verticali, sinonimo di eccessiva compressione. La facciata era soggetta a 
meccanismo di ribaltamento fuori dal piano e la cupola interna in gran parte collassò. 
Dopo il sisma è stato installato un presidio temporaneo al fine di prevenire ulteriori crolli. 


La chiesa è stata restaurata e riaperta al pubblico nel 2018. 
Descrizione della sperimentazione effettuata 


L'attività di monitoraggio è iniziata nel novembre 2009 e terminata nel febbraio 2013 da 


parte dell’Università di Venezia. 


In Tabella 6 si riportano le caratteristiche del sistema di monitoraggio. 
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Tabella 6: Caratteristiche principali del sistema di monitoraggio dinamico. 


19 ACCELEROMETRI PIEZOELETTRICI PCB 393C MONODIREZIONALI [AM] 


3 ACCELEROMETRI PCB 356B18 TRIDIREZIONALI [AT] 
24 schede di acquisizione DYMAS24-ALBEN 
2 VELOCIMETRI Lennartz Le3D-s [AT-Ingv] 
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Figura 12: Disposizione planimetrica della strumentazione adottata. 
In Figura 12 viene illustrata la disposizione planimetrica della strumentazione adottata. 


Inizialmente sono stati adottati 19 accelerometri monodimensionali (AM) e 3 
accelerometri tridimensionali (AT). Durante l'ultimo mese di sperimentazione (dal 18 
gennaio al 13 febbraio 2013) sono stati aggiunti due velocimetri (AT-Ingv). Questi 
strumenti sono notevolmente sensibili e sono capaci di rilevare frequenze di risonanza 


molto basse. 
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È stata calcolata l'ampiezza degli spettri di Fourier (FAS Fourier Amplitude Spectra) per le 
frequenze acquisite dagli strumenti. Gli spettri delle accelerazioni elaborati mostrano una 


tendenza decrescente con il crescere della frequenza. 


Durante il periodo di installazione del monitoraggio, sono stati registrati diversi eventi 
sismici tra cui due terremoti locali di magnitudo 3.6Mw a 11km a est della città e 4.8Mw 


a 70kma sud. Peri due eventi sono state calcolate: 
e L’accelerazione di picco (PA); 
e La velocità (PV); 
e Lospostamento (PD). 


Risultati e conclusioni 


Da questo studio è stato possibile evidenziare il confronto tra gli accelerometri e i 
velocimetri. Difatti, come mostrato in Figura 13 si può notare che gli accelerometri non 
colgono frequenze inferiori a 2 Hz dove la componente di rumore è predominante. | 
velocimetri invece sono molto più sensibili nel campo delle basse frequenze e arrivano ad 
indagare fino a valori di frequenza di 0,2 Hz. In dettaglio, questi strumenti riescono a 
cogliere la frequenza di risonanza di 0,5 Hz alla base della chiesa che corrisponde con la 
frequenza di risonanza del sottosuolo del centro dell’Aquila. (Di Giulio, Vassallo, Boscato, 


Dal Cin, & Russo, 2014) 
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Figura 13: Comparazione del rumore registrato tra velocimetri (Nero) e Accelerometri (Rosso). 
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Dalla sperimentazione effettuata, si è potuto notare inoltre che se questi strumenti 
vengono disposti a terra, la misura può essere affetta da un’elevata componente di 


rumore. 


Dall'analisi condotta sugli eventi sismici registrati, emerge che gli spostamenti, le velocità 
e le accelerazioni mostrano una tendenza crescente con l'altezza. Questa tendenza è 


chiaramente dimostrata nei grafici riportati in Figura 14 relativi al sisma di magnitudo 4.8. 
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Figura 14 Terremoto di magnitudo 4.8. Nei grafici di destra vengono mostrati gli spostamenti, le velocità e le 
accelerazioni di alcune stazioni di misura al crescere dell'altezza, mentre nei grafici di sinistra vengono plottate 
le stazioni accelerometriche alle diverse altezze indicando con il rombo la direzione x con il triangolo la 
direzione y e con il cerchio la direzione z. 

Per rappresentare le frequenze di risonanza della struttura è stata usata l’analisi tempo- 
frequenza che risulta essere affidabile quando la struttura presenta un comportamento 


non stazionario e variabile neltempo come appunto quello indotto da un terremoto. 


In conclusione, è stata calcolata anche la velocità delle onde di taglio che è pari a 150- 


200m/s. 
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1.2.1.4. La Chiesa di San Pietro a Coppito in L'Aquila 


Figura 15: Chiesa di San Pietro a Coppito in L’Aquila. 


Cenni storici 


La Chiesa venne edificata nella metà del XIII secolo probabilmente sul luogo di un 


preesistente tempio. 


Il monumento venne distrutto una prima volta in seguito al terremoto dell'Aquila nel 1703 
dove crollarono sia la facciata che la torre campanaria. La sua ricostruzione richiese 
moltissimi anni. In seguito agli eventi sismici che hanno colpito l'Aquila nel 2009, la chiesa 
di S. Pietro è stata di nuovo severamente danneggiata. Principalmente, si sono innescate 
le seguenti fenomenologie: la prima ha interessato la facciata che ruotando verso 
l'esterno ha causato il collasso della parte sommitale destra; peraltro, il sisma ha causato 


sia il crollo del campanile che di una parte dell'abside. 
I lavori di ricostruzione sono ancora in corso. 


Descrizione della sperimentazione effettuata 


In via preliminare, è stata condotta un'intensa campagna d’indagine in sito, riportata in 


Figura 16 che ha previsto la realizzazione delle seguenti prove: 
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e Esecuzione di test sonici diretti, indiretti e semi diretti; 


e Prelievo di campioni in alcune porzioni della struttura 
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Figura 16: Piano di indagine. 


Inoltre, sono state effettuate delle prove con i martinetti piatti. Un primo test è stato fatto 
a marzo 2011 e un secondo test a settembre 2011. In particolare, i test della prima fase 
hanno avuto l’obiettivo di caratterizzare la muratura senza portarla allo sforzo ultimo, 
mentre i test effettuati nella seconda fase sono serviti per valutare gli effetti causati dagli 


effetti ambientali. 


È stato eseguito anche il monitoraggio dinamico attraverso l’uso di 12 accelerometri. 
Questi sono stati posizionati alle altezze di 6 e 10 metri. In Figura 17 si riporta nel dettaglio 


la disposizione planimetrica della rete di sensori adottata. 


Figura 17: Disposizione del sistema di monitoraggio dinamico. 
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Per l'identificazione dei parametri modali è stato usato il software ARTeMIS. 


Risultati e conclusioni 


| test non distruttivi condotti sugli elementi della navata centrale e dell’arco trionfale 
rivelano un differente stato di compattezza della sezione che quindi si traduce in un 
comportamento asimmetrico degli archi. Attraverso i test sonici è stato possibile 
identificare un danno diffuso in entrambe le direzioni, le quali hanno confermato la 
presenza di zone vuote e caratterizzate da differenti proprietà elasto-plastiche tra la zona 
laterale e quella centrale. Sui campioni prelevati sono state condotte delle prove che 
mostrano una buona resistenza a compressione (Figura 18) per quanto concerne i pilastri 
interni e la torre campanaria, viceversa, la prova condotta sui provini della facciata 


esterna denota una scarsa qualità dei materiali. 
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Figura 18: Diagramma sforzo-deformazione dei campioni prelevati nella chiesa. 


Inoltre, in Tabella 7, si riportano i risultati delle prove dei martinetti piatti. 
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Tabella 7: Risultati test Martinetti piatti. 
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Le prove di vibrazione ambientale condotte sulla chiesa seppur danneggiata hanno 
permesso di identificare il comportamento dei principali macroelementi e l'interazione 
tra i principali sistemi strutturali. Emerge che l’articolata configurazione plano-altimetrica 
amplifica gli effetti delle azioni dinamiche. La snellezza della struttura e l'assenza di 
elementi di irrigidimento ha reso l’edificio molto deformabile e soggetto quindi a 
ribaltamenti fuori dal piano e meccanismi di traslazione che giustificano in parte lo stato 


di danno emerso in conseguenza del sisma. (Boscato, Dal Cin, Russo, & Sciarretta, 2014) 


In Figura 19 sono riportate le frequenze relative ai primi dieci modi della struttura. Si nota 
che a causa delle irregolarità in pianta e in altezza e della diversa distribuzione delle 
rigidezze e della geometria della struttura il primo modo è caratterizzato da una 
deformata modale di tipo torsionale, mentre la torre campanaria nonostante fosse 


collegata alla chiesa non è coinvolta nella risposta strutturale globale della struttura. 
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Figura 19: Forme modali della chiesa. 


Emerge quindi che il comportamento globale della chiesa è governato dalla navata, dalla 


facciata, dall'arco trionfale e dalla complessa articolazione del doppio transetto e dei tre 


absidi. La torre campanaria risulta il terzo elemento del sistema per ragioni legate a 


rigidezza della struttura, condizioni al contorno e alla geometria. 
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1.2.1.5. La Basilica di San Nicola da Tolentino (Mc) 


Figura 20: Basilica di San Nicola da Tolentino. 


Cenni storici 


La Basilica di San Nicola da Tolentino inizialmente intitolata a San Giorgio venne costruita 


tra il XIII e il XIV secolo per poi essere completata nei secoli successivi. 


Gli eventi sismici susseguitesi nel 2016 hanno rilevato danni distribuiti su tutta la struttura 
della Basilica. Significative sono le lesioni in corrispondenza dei contrafforti della facciata 
a livello della copertura. Il movimento di distacco della facciata è evidenziato anche nella 
parte inferiore da fessure in corrispondenza della cornice decorativa interna su entrambi 
i lati della navata. In seguito all'entità delle lesioni, è stato installato un sistema di 
monitoraggio per evitare la caduta di elementi murari e decorativi. In corrispondenza 
dell'abside sono presenti delle lesioni che indicano una probabile lesione del corpo 
centrale. Fessurazioni diffuse sono presenti anche sugli archi laterali della navata, sulla 


cappella confessionale e sul cappellone. 
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Descrizione della sperimentazione effettuata 


Nella struttura è stato predisposto un sistema di monitoraggio con lo scopo di acquisire il 
comportamento dinamico globale dei macroelementi principali (facciata, navata, abside 
e torre campanaria) danneggiati durante la sequenza sismica in modo da valutare 


l'attivazione o l'evoluzione dei meccanismi di danno individuati. 


Le caratteristiche del sistema di monitoraggio sono riportate in Tabella 8, e la disposizione 


della strumentazione è evidenziata in Figura 21. 


Tabella 8: Caratteristiche principali del sistema di monitoraggio dinamico. 


20 ACCELEROMETRI MONOASSIALI 
3 TERMOCOPPIE 
SISTEMA DI ACQUISIZIONE 24 BIT 


Figura 21: Posizione dei sensori acceleromettrici (in giallo) e delle termocoppie (in rosso). 


La disposizione degli accelerometri è basata sul rilievo dello stato di danno. Pertanto, sono 
state collocate due coppie di sensori sulla facciata e sui muri longitudinali della chiesa per 
controllare l’efficacia del collegamento e il loro comportamento sotto l’effetto di 


eventuali eventi sismici. (Sabia, Aoki, Costanzo, Lancellotta, & Quattrone, 2017) 


24 


Le time history sono state acquisite in intervalli temporali di un’ora per un totale di nove 
settimane. Sulla base di questi segnali sono stati poi identificati i parametri modali 
attraverso l’algoritmo SSI. | modi reali della struttura sono stati individuati applicando i 


seguenti criteri di stabilizzazione: 
e Variazione della frequenza tra ordini successivi inferiore all'1%; 
e Variazione dello smorzamento tra ordini successivi inferiore al 10%; 
e Esclusione dei modi con € > 10% e € < 0%; 
e MAC>95%. 


È stato realizzato un monitoraggio continuo e la rappresentazione nel dominio tempo- 
frequenza che ha permesso di valutare l'evoluzione nel tempo delle frequenze naturali 


dei principali modi di vibrare. 
Risultati e conclusioni 


I dati acquisiti nei primi due mesi del monitoraggio hanno permesso l’identificazione dei 


modi di vibrare della struttura associati ai macroelementi maggiormente danneggiati. 


In Figura 22 si riportano i diagrammi di stabilizzazione e di clusterizzazione dei risultati al 


variare dell’ordine del sistema. 
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Figura 22: Diagrammi di stabilizzazione e clustering. 
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I primi due modi fanno riferimento alla torre campanaria, il terzo ed il quarto sono i modi 
globali della basilica nelle due direzioni orizzontali, mentre il quinto individua il modo di 


vibrare della facciata. 


Attraverso l'installazione di un sistema di monitoraggio permanente in futuro si potranno 
valutare eventuali deviazioni del comportamento strutturale in base ai parametri acquisiti 
in questa prima fase. Inoltre, potrà essere valutata la risposta della struttura durante 


eventuali eventi sismici. 


1.2.2. Prove dinamiche applicate a edifici di pregio 


In questo paragrafo, vengono riportati i casi studio in cui sono state realizzate delle prove 
dinamiche su edifici di importante valenza storico-architettonica. Gli aspetti indagati sono 


i medesimi della sezione precedente. 


1.2.2.1. Il Duomo di Milano 


Figura 23: Il Duomo di Milano. 
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Cenni Storici 


II Duomo di Milano è noto in tutto il mondo per essere uno dei più grandi monumenti ad 


esser stato costruito e per la sua unica architettura. 


L'inizio della costruzione avvenne nel 1386, fu promossa dal vescovo Antonio de' Saluzzi 
e quindi autorizzata dal duca di Milano Gian Galeazzo Visconti. Per la realizzazione della 
chiesa vennero fatte abbattere prima la Cattedrale di Santa Maria Maggiore e in un 


secondo momento la Basilica di Santa Tecla. L’edificazione venne conclusa nel 1813. 


La cattedrale presenta pianta a croce latina con dimensione 66x158m. Negli ultimi anni 
del XX secolo al fine di ridurre la spinta sulle pareti ortogonali sono stati installati 122 


tiranti d'acciaio all'altezza dei capitelli dei pilastri. 


Descrizione della sperimentazione effettuata 


Nella Cattedrale di Milano è stato predisposto un sistema di monitoraggio statico e un 


sistema di monitoraggio dinamico. 


Il sistema di monitoraggio statico è completamente computerizzato ed è consistito nelle 


misurazioni delle seguenti quantità: 
e  Sforzi quasi statici attraverso l’utilizzo di estensimetri wireless a cavo vibrante; 


e Rotazioni quasi statiche misurate in cima ad alcuni pilastri selezionati 


arbitrariamente e su tre livelli nella guglia principale; 
e Parametri ambientali interni ed esterni; 
e  Velocitàsutre livelli nella guglia principale. 


In Tabella 9 vengono riportate le caratteristiche del sistema di monitoraggio statico e in 


Figura 24 si riporta la disposizione plano-altimetrica dei sensori adottati. 
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Tabella 9: Caratteristiche della strumentazione adottata per il monitoraggio statico. 


12 ESTENSIMETRI A FILO VIBRANTE 


em 


Figura 24 : Mappa dei dispositivi installati all’interno della Cattedrale di Milano e sulla guglia principale. 


Per quanto concerne il monitoraggio dinamico, sono stati utilizzati dei sismometri cioè dei 


trasduttori elettrodinamici di velocità e la loro scelta è stata motivata da: 
e Alta sensibilità ed eccellenti performance nel campo delle basse frequenze; 
e Alimentazione non necessaria; 


e Permettono di ottenere una buona stima degli spostamenti dall’integrazione 


delle velocità. 
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In Tabella 6 si riportano le caratteristiche del sistema di monitoraggio dinamico e in Figura 


25 si riporta la disposizione plano-altimetrica dei sensori adottati. 


Tabella 10: Caratteristiche del sistema di monitoraggio dinamico. 


13 SISMOMETRI BIASSIALI (sopra pilastri) 
1 SISMOMETRO MONOASSIALE (sopra pilastri) 
9 SISMOMETRI MONOASSIALI (guglia principale) 
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Figura 25: Mappa dei dispositivi installati all’interno della Cattedrale di Milano. 


I parametri modali della cattedrale e della guglia principale sono stati estratti 
indipendentemente e sono stati processati mediante l'algoritmo SSI-COV. È stata adottata 
una procedura di identificazione automatica dei parametri modali che viene effettuata in 


due fasi: 


e La ModalParameter Estimation (MPE) è stata condotta attraverso una procedura 
automatica di interpretazione dei diagrammi di stabilizzazione basata sulla 
sensibilità delle frequenze e delle forme modali al cambiare dell'ordine del 


modello. Per ogni set di dati sono stati eliminati i modi spuri dal controllo del 
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rapporto di smorzamento e dall’indice di complessità modale (MCI). | poli che 
hanno mostrato frequenze e forme modali simili sono stati clusterizzati 
determinando per ogni cluster i parametri modali (frequenze di risonanza, forme 


modali e rapporti di smorzamento). (Gentile, Ruccolo, & Canali., 2019) 


e Il Modal Tracking (MT) è stato condotto al fine di valutare l'evoluzione dei 
parametri modali di ogni modo ed è basato sulla variazione delle frequenze e della 


MAC rispetto ad una lista di modi di base. (Gentile, Ruccolo, & Canali., 2019) 


Dall'ottobre del 2018 a febbraio del 2019 è stato condotto un monitoraggio dinamico. 
Durante questo arco temporale è stata fatta un’acquisizione all’ora e i dati sono stati 
processati automaticamente secondo la procedura descritta in precedenza. Si è potuta 
valutare la correlazione tra le frequenze della struttura e i parametri ambientali come 


temperatura e umidità. 
Risultati e conclusioni 


In seguito all’installazione del sistema di SHM applicato al Duomo di Milano è stato 
possibile innanzitutto constatare l‘efficacia di strumenti sensibili per l’analisi modale 
operativa permettendo un'accurata stima dei parametri modali della struttura, in 
particolare ciò ha consentito l’individuazione di otto frequenze di risonanza tra 0 e 5 Hz 


come mostrato in Figura 26. 


30 


idrata drenante nenti fos 
4.0 
— 3.5 
N 
E for 
> 3.0 
2 fcs 
: fa 
© 25 fc 
© 
dl 
i 


10 
16/10/2018 15/11/2018 15/12/2018 14/01/2019 13/02/2019 
measurement date 


Figura 26: Frequenze di risonanza individuate. 


Anche se il periodo di osservazione è stato limitato a quattro mesi, la correlazione 
frequenza-temperatura rivela una tendenza distintiva che va in antitesi con molti studi a 
lungo termine su strutture in muratura. Come si può notare in Figura 27 tutte le frequenze 
di risonanza della Cattedrale di Milano crescono più o meno linearmente con il decrescere 
della temperatura e la tendenza negativa delle frequenze modali sulla temperatura è 


plausibilmente determinata dagli effetti attivi esercitati dai tiranti metallici. 


Figura 27: Correlazione tra temperatura esterna e le prime quattro frequenze misurate. 
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Infine, si è scoperto che le forme modali della cattedrale non esibiscono apprezzabili 
fluttuazioni con gli effetti delle vibrazioni ambientali e questo potrebbe portare a 


considerare nei successivi studi di SHM questo parametro come invariante. 


1.2.2.2. Il Santuario di Vicoforte a Cuneo 


Figura 28: Santuario di Vicoforte a Cuneo. 


Cenni storici 


Il Santuario di Vicoforte è una struttura in muratura la cui costruzione ebbe inizio nel 1596 
ed è proseguita a fasi alterne fino ad essere completata nel 1884. L’edificio si 
contraddistingue per la sua cupola a sezione orizzontale ellittica che risulta essere la più 


grande di tale forma al mondo. 


Il tamburo della cupola ha sofferto nel corso degli anni di significativi problemi strutturali 
causati in parte dall'aggiunta di masse, che quindi hanno causato anche dei cedimenti nel 
terreno, ma in larga parte anche a causa della configurazione strutturale del duomo 
stesso. Nel 1983 a causa del peggioramento del quadro fessurativo dopo un'accurata 
campagna d'indagine sono stati installati 56 tiranti attivi che vennero posizionati 


all'interno di fori praticati nella muratura nella parte superiore del tamburo, lungo 14 
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tangenti intorno al perimetro a 4 livelli differenti come mostrato in Figura 29. Nel 1997 i 
tiranti sono stati ritesati per compensare le perdite da stress. (Chiorino, Spadafora, 


Calderini, & Lagomarsino, 2008) 


a} Wost 


Figura 29: a) disposizione planimetrica dei tiranti all’interno del tamburo; b) dettaglio dei tiranti. 


Descrizione della sperimentazione effettuata 


In seguito all'installazione dei tiranti metallici è stato posto in opera un sistema di 
monitoraggio statico per valutare gli effetti del sistema di rinforzo e la propagazione delle 


fessure. 


Successivamente, nel dicembre del 2015, anche un sistema di monitoraggio dinamico 
permanente è stato installato sul sistema tamburo-cupola-lanterna le cui caratteristiche 


vengono sintetizzate in Tabella 11 e la disposizione plano altimetrica in Figura 30. 


Tabella 11: Strumentazione adottata per il sistema di monitoraggio dinamico. 


12 ACCELEROMETRI PIEZOELETTRICI (PCB 393B12) 
UNITA' DI ACQUISIZIONE 
UNITA' MASTER 
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Figura 30: Disposizione degli accelerometri nelle sezioni A-A e B-B. 


Il sistema di acquisizione è stato progettato con uno schema di acquisizione master/slave 
per limitare la distorsione causata dalla lunghezza dei cavi sopra ai 50 m di lunghezza e 
per ridurre l’invasività. | dati acquisiti dagli accelerometri nella cripta sono trasmessi 
all'unità slave e quindi all'unità master. | parametri registrati dagli altri strumenti 
convergono direttamente all'unità master e il tempo di acquisizione viene poi 


sincronizzato con il sistema GPS. 


Per l'elaborazione dei dati del sistema di monitoraggio è stata adottata la tecnica della 
cointegrazione, che si basa sugli algoritmi SVM e RVM, al fine di rimuovere i dati non 


stabili. 


È stata effettuata anche un’analisi dei parametri ambientali per cercare di comprendere 


eventuali variazioni e correlazioni con le frequenze misurate della struttura. 


AI fine di validare l'efficacia del metodo di rilevamento del danno è stato simulato uno 


scenario di danneggiamento su un modello agli elementi finiti del santuario che era 
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disponibile da precedenti ricerche. Il comportamento meccanico dei materiali è stato 
assunto elastico lineare. Nelle zone con maggiori sollecitazioni, riportate in arancio in 


Figura 31, per simulare lo stato di danno è stato ridotto il modulo di Young del 40%. 


(a) (b) 
Figura 31: Modello del Santuario di Vicoforte nelle due viste Sud Nord ed Est Ovest. 
Risultati e conclusioni 


Il sistema di monitoraggio ha permesso l'individuazione di sedici frequenze che vengono 
riportate in Figura 32. Si può notare come le prime due frequenze siano molto più regolari 


rispetto alle altre. 
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Figura 32: Frequenze identificate del santuario. 
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Lo studio condotto mostra i vantaggi dell’uso della tecnica della cointegrazione che è 
basata su un’idea semplice e versatile e quindi può essere usata con diversi sistemi con 
schemi statici, materiali ed età differenti e può essere anche implementata senza alcun 
dato sulle condizioni ambientali. La cointegrazione applicata alle frequenze identificate 


dimostra che non ci sono danni nel periodo indagato. (Coletta, et al., 2018) 


In Figura 33 sono riportati tutti i parametri ambientali misurati nel periodo di 


funzionamento del sistema di monitoraggio. 
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Figura 33: Prime due frequenze e parametri ambientali registrati nei 4 mesi di monitoraggio. 


In particolare, si scorge una tendenza anomala attorno alla 170° osservazione, dove si 
hanno condizioni di neve e temperatura inferiore allo 0 e un temporaneo innalzamento 
delle prime due frequenze della struttura. Probabilmente in questo caso l’azione del 


ghiaccio ha prodotto un effetto di irrigidimento della chiesa. 


Un’ultima osservazione è stata effettuata sulla forma della cupola che essendo ellittica 
porta ad avere due frequenze diverse e molto ravvicinate tra loro. Nel caso in cui la forma 


fosse stata circolare le prime due frequenze rilevate sarebbero dovute coincidere. 
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1.2.2.3. Il Battistero della Chiesa di San Giovanni a Firenze 


TOLD 


a 


Figura 34: Il battistero della Chiesa di San Giovanni a Firenze. 
Cenni Storici 


Il Battistero di San Giovanni a Firenze è considerato un eccezionale monumento a livello 
mondiale che venne consacrato nel 1059 e sorge di fronte alla Cattedrale di Santa Maria 


del Fiore. La struttura presenta una tipica pianta ottagonale con un diametro di 25,60 m. 


Dalla sua costruzione il monumento ha subito diversi sismi di varia magnitudo (1-VIII MCS 
scale). Negli ultimi anni ha subito anche l'influenza del traffico veicolare, che dall'ottobre 


2009 ha portato ad una chiusura della zona adiacente. 


Descrizione della sperimentazione effettuata 


È stato effettuato un monitoraggio dinamico per 42 ore dal 14 al 16 ottobre 2013 con 10 
stazioni sismiche le cui caratteristiche sono riportate in Tabella 12. Ogni stazione sismica 


è stata dotata di sismometri a tre componenti. 
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Tabella 12: Caratteristiche del monitoraggio dinamico. 


6 SISMOMETRI LENNARTZ 3D/5 s (sensibilità 400 V/m/s funzione di trasferimento 5s) 


3 SISMOMETRI GURALP CMG-ET (sensibilità 2400 V/m/s con funzione di trasferimento 105) 
1 SISMOMETRO GURALP CMG-40T (sensibilità 800 V/m/s con funzione di trasferimento 305) 
SCHEDE CONVERSIONE SEGNALE A/D 24 BIT (GURALP CMG24) 


In Figura 35 viene mostrata la disposizione delle stazioni sismiche, nel dettaglio: 


e la stazione BO1 è stata collocata al piano terra e usata come stazione di 


riferimento al suolo; 


e lestazioni B02-B03-B04-BOS5 sono state posizionate alla base del Matroneo; 


e lestazioni B06-B07-B08-B09 sono state posizionate al di sopra del Matroneo; 


e LaB10, si trova sul lucernario. 


10m N 


# B10 - (Lennartz 55) 


First Women's Gallery Second Women's Gallery 


Figura 35: Disposizione delle 10 stazioni sismiche nel Battistero di San Giovanni. 
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In questo studio i dati derivanti dai segnali sono stati processati usando il metodo EFDD 
che è un'estensione del Peak picking ed è basato sulla definizione dei valori singolari e 
della Power Spectral Density Matrix. Inoltre, è stata usata una procedura di identificazione 


automatica introdotta da Brinker, per analizzare i parametri modali della struttura. 


In conclusione, è stata realizzata una modellazione 3D per confrontare i risultati con i dati 
sperimentali. Inoltre, è stata fatta una valutazione dello spostamento corrispondente ad 
ogni forma modale attraverso un’integrazione della velocità sismica al suolo e un filtraggio 


a 1Hz attorno alla frequenza naturale di ogni modo. 


Risultati e conclusioni 


L'applicazione della procedura automatica basata sull’algoritmo EFDD per il calcolo dei 
parametri modali si è rilevata efficace portando alla definizione di sei forme modali 
fondamentali che ampliano la conoscenza di questo monumento antico. Il metodo si è 
rivelato efficace anche perché riesce a cogliere bene i primi due modi che sono separati 
solo da 0,21 Hz l’uno dall’altro. Inoltre, come mostrato anche in Figura 36 a), dall’analisi 
dei dati delle 42 ore di monitoraggio è emerso che i primi cinque modi sono sempre 


presenti, mentre il sesto sembra esserci solo nel periodo diurno e questo potrebbe essere 


dovuto all'attività antropica nelle zone circostanti. 


Si può osservare sempre dalla stessa figura, una leggera crescita dei valori delle frequenze 
attorno alle ore 12:00, che può essere relativa ad un incremento degli effetti della 
temperatura. Questa correlazione è stata studiata in ricerche a lungo termine che hanno 
mostrato come la chiusura delle fessure sia causata dall’espansione termica con 
conseguente crescita della rigidezza e delle frequenze modali. Per confermare tale 
tendenza sono necessari ulteriori monitoraggi a lungo termine. (Lacanna, Ripepe, 


Marchetti, Coli, & Garzonio, 2016) 
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Damping Ratio (%) 


splacement (um) 


Di 
T=] 
o 
9 


00.00 06:00 12.00 18:00 00:00 
2013/10/15 time 2013/10/16 


Figura 36: a) Registrazioni tempo-frequenze; b) Registrazioni tempo-rapporti di smorzamento; c) 
Registrazioni tempo-spostamenti. 


Gli spostamenti, in Figura 36c), variano da un minimo di 0,08 um di notte a un massimo 
di 0,7 um di giorno, sinonimo della presenza consistente dell’attività antropica nella zona 


del Battistero. 


In Tabella 13 sono mostrate le caratteristiche dinamiche del battistero in termini di 


frequenze, rapporti di smorzamento e spostamenti. 
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Tabella 13: Risultati di frequenze smorzamenti e spostamenti relativi alle forme modali. 


Mode number Frequency (Hz) 
Mean Standard deviation 

1 2.67 0.010 

2 2.88 0.017 

3 4.80 0.029 

4 5.88 0.020 

5 6.34 0.029 

6 693 0.061 
Damping ratio (*) Displacement (pm) 
Mean Standard deviation Mean 
141 0.166 0.247 
141 0.200 0.221 
1.64 0.306 0.071 
1.23 0.170 0.049 
1.56 0.220 0.048 
1.05 0.210 0.035 

1.2.2.4. La Cattedrale di Colonia in Germania 


Figura 37: La Cattedrale di Colonia. 
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Cenni storici 


La Cattedrale di Colonia è un grandioso tempio in stile gotico la cui costruzione è durata 
più di 600 anni, essendo iniziata nel 1248 e terminata solo nel 1880. La cattedrale venne 
fortemente colpita dai bombardamenti della Seconda guerra mondiale, ma 
fortunatamente non collassò ed i lavori di restauro terminarono nel 1956. Nel 1996 è stata 


inserita nell'elenco dei patrimoni dell'umanità dell'Unesco. 


La struttura in esame è soggetta a notevoli vibrazioni ambientali, infatti vicino ad essa è 
presente la stazione centrale nella quale ogni giorno transitano circa 1200 treni, una 
strada molto trafficata passa a Nord dell'edificio e un mega parcheggio si trova a Sud della 
stessa. A causa di ciò le vibrazioni ambientali cui è soggetto l'edificio sono di gran lunga 


superiori durante il giorno piuttosto che durante la notte. 


Descrizione della sperimentazione effettuata 


L'obiettivo principale del controllo della Cattedrale di Colonia è quello di monitorare il 
comportamento dinamico durante i terremoti locali. In aggiunta, il monitoraggio continuo 


è anche utile per valutare l’attività antropica e l’effetto delle condizioni climatiche. 


A tal proposito, sono state installate cinque stazioni sismiche le cui caratteristiche sono 


riportate in Tabella 14 e la disposizione plano-altimetrica in Figura 38. 


Quattro di queste stazioni hanno sensori a tre componenti, mentre la stazione BA21A è 


dotata solo di due componenti in direzione orizzontale. 


Tabella 14: Caratteristiche principali del sistema di monitoraggio dinamico. 


5 SISMOMETRI A 2-3 COMPONENTI KINEMETRICS 


SCHEDE CONVERSIONE SEGNALE A/D A 4 CANALI 
NANO SERVER PC 
STAZIONE METEO (nelle vicinanze della struttura) 
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Figura 38: Disposizione delle stazioni sismiche nella cattedrale, all’interno del cerchio è indicato il numero di 
componenti di ogni stazione. 


Sono state valutate le frequenze proprie della struttura su un periodo di sette giorni (in 
assenza di vento e condizioni di alta pressione) e ogni quattro ore sono state analizzate 
finestre temporali di quattro minuti che hanno fornito un range di frequenze tra 0,2 Hz e 


2,0 Hz. 
Con i dati metereologici acquisiti nelle vicinanze della struttura è stato possibile studiare: 
e Lacorrelazione tempo-frequenza per tutto il 2010; 


e Lacorrelazione frequenza-velocità del vento nel periodo da gennaio a marzo 2010 
attraverso una stazione anemometrica situata a 100m da terra. Nella torre nord 
si sono registrate velocità medie del vento comprese in un range da 2,5 a 47,5 


m/s; 


e La correlazione frequenza-temperatura. 
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È stata oggetto di studio anche la sollecitazione dovuta al suono delle campane, infatti 
nella cattedrale è presente una gigantesca campana, chiamata dagli abitanti di Colonia 


Petersgloke, dal peso di 24000 kg con un batacchio di 800 kg. 


Risultati e conclusioni 


Il sistema di monitoraggio della struttura ha portato a rilevare quattro frequenze proprie. 
Particolari conclusioni possono essere tratte in merito alle correlazioni citate in 


precedenza. 


Di particolare interesse è stata la correlazione frequenza — velocità del vento. In dettaglio, 
si nota che l'accelerazione spettrale relativa ai primi tre modi cresce esponenzialmente 
con la velocità del vento come si evince dal grafico di Figura 39. Si nota che le frequenze 
proprie sono meno variabili della velocità del vento, che può esibire cambiamenti 
transitori e quindi possono permettere di definire il livello di allerta superato il quale può 
determinarsi la chiusura delle aree circostanti e della chiesa per evitare il rischio del crollo 


di alcune parti. (Hinzen, Fleischer, Schock-Werner, & Schweppe, 2012) 
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Figura 39: Andamento della velocità media del vento con l'accelerazione spettrale. 


Per quanto concerne la correlazione frequenza temperatura, è stato osservato un 


cambiamento nei mesi invernali soprattutto della frequenza relativa al quarto modo che 
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è aumentata passando da 1,4 Hz a 1,5 Hz (Figura 40). Sempre in merito a questa 
correlazione, in generale si è osservato un aumento delle frequenze con la diminuzione 
della temperatura dovuto probabilmente all'umidità che occupando i pori dell’arenaria, 


gelando, causa un incremento della rigidezza in grado di far cambiare la frequenza. 
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Figura 40: Correlazione frequenza temperatura. 


Anche gli effetti del suono della grande campana, ha condotto a riflessioni interessanti. A 
causa dell’elevata massa il movimento sinusoidale della campana è in grado di produrre 
vibrazioni misurabili in tutte le cinque stazioni della cattedrale. In Figura 41, sono riportati 
gli spostamenti registrati in una finestra temporale di due minuti da quando la campana 
inizia a suonare. Si nota che le frequenze principali di questi movimenti sono in accordo 


con la prima e la seconda frequenza per le componenti NS e EW. (Hinzen, Fleischer, 


Schock-Werner, & Schweppe, 2012) 
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Figura 41: Spostamenti registrati nelle stazioni BA19 BA20 BA21 quando la campana inizia a suonare. 


1.2.2.5. La Chiesa di San Torcato a Guimaràes in Portogallo 


Figura 42: Chiesa di San Torcato a Guimaràes. 
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Cenni storici 


La Chiesa di San Torcato si trova nella città di Guimaràes nel nord del Portogallo. La 
costruzione dell’edificio ebbe inizio nel 1871 e venne completata nel 2008. A causa del 
lungo processo di edificazione sono stati usati diversi materiali. Le torri e i muri della 
navata principale sono costituiti principalmente in muratura con blocchi di granito, 
mentre l'abside e l’altare principale sono stati realizzati in cemento armato e rivestiti di 
granito. Negli ultimi anni della sua costruzione la chiesa è stata oggetto di monitoraggio 
in quanto nella facciata, soprattutto nella chiave di volta dell'arco principale e sulle due 
torri era presente una fessurazione diffusa con ampiezza di circa 5 cm sul fronte e al livello 


del coro erano superiori a 3cm. 


Descrizione della sperimentazione effettuata 


Nel 1998 è stato installato nella chiesa un semplice sistema di monitoraggio low-cost. La 


strumentazione adottata è riportata in Tabella 15. 


Tabella 15: Caratteristiche principali del sistema di monitoraggio statico. 


FESSURIMETRI 


SET DI 3 PIASTRE D'ACCIAIO (Riferimento per le misure di spostamento) 


STAZIONE TOTALE (perle misura dell'inclinazione e degli spostamenti delle torri) 


MISURA DELLA TEMPERATURA SULLA SUPERFICIE DI OGNI PUNTO OGGETTO DI STUDIO 


Dai risultati si nota che il comportamento delle fessure è fortemente influenzato dalle 


condizioni ambientali. 


Gli spostamenti della torre sono stati valutati usando una stazione di misura totale 
prendendo le misure lungo l'altezza, mentre le misure inclinometriche sono state fatte 
con un piombo posizionato all’interno di una torre. Uno schema delle indagini viene 


riportato in Figura 43. 
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Figura 43: Dettaglio delle indagini statiche. 


È stato realizzato anche un rilievo geotecnico basato sulla realizzazione di 31 CPT e la 
realizzazione di due sondaggi aperti per osservare direttamente le fondazioni e prelevare 


i campioni indisturbati. 


Nel maggio 2009 è stata predisposta un’ampia campagna di identificazione dinamica 
globale basata sull'analisi modale operativa OMA. | sensori di riferimento sono stati 
posizionati in 35 punti riportati in Figura 44. A causa del numero di accelerometri usati e 
delle lunghe distanze, le misurazioni sono state effettuate in nove configurazioni di prova 


utilizzando quattro punti di riferimento. (F.Ramos, Aguilar, Lourenco, & Moreira, 2013) 


VLLLLLIILILLIZZIIZZI 


Figura 44: a) Layout dei sensori durante l’identificazione modale; b) Sistema di acquisizione sopra la navata; 
c) Acceleromettri. 
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| parametri modali della struttura sono stati ottenuti usando gli algoritmi FDD, (Frequency 
Domain Decomposition) EFDD (Enhanced Frequency Domain Decomposition) e CC-SSI 


(Cristal Clear Stochastics Sub-Space Identification). 


È stato realizzato un modello agli elementi finiti della basilica. A causa della complessità è 
stata modellata solo la facciata e una parte della navata centrale. | risultati del modello 
analitico sono sensibili delle proprietà dei materiali. Per effettuare la procedura di 
updating sono stati usati i risultati derivanti dall’identificazione modale e delle indagini 


svolte preliminarmente. 


Nella chiesa è stato installato nel 2009 un sistema di monitoraggio dinamico, le cui 
caratteristiche sono riportate in Tabella 16, per cercare di capire se ci fosse un 


cambiamento della risposta dinamica globale nel tempo. 


Tabella 16: Caratteristiche principali del sistema di monitoraggio dinamico. 


4 ACCELEROMETRI PIEZOELETTRICI 


DAQ (Risoluzione 24 bit con 4 canali di misura) 


UPS (Uninterruptible Power Supply) 


COMPUTER 
MODEM 


Per analizzare la grande mole di dati è stato usato l'algoritmo SSI Data-Driven. 


L’identificazione modale è condotta a partire dalla costruzione dei diagrammi di 
stabilizzazione e dalla loro elaborazione. In particolare, per eliminare i trend spuri e 
trovare quindi le frequenze proprie della struttura è stata adottata una procedura di 
clustering non gerarchica che usa l'algoritmo K-means. Tuttavia, questa procedura 
dipende fortemente dalla selezione del cluster di dati iniziale. L'algoritmo proposto 
considera due step consecutivi: la costruzione del diagramma di stabilizzazione e l'ordine 
del modello per il quale i parametri dovranno essere stimati. L'obiettivo di ogni livello del 
cluster è quello di considerare frequenze e rapporti di smorzamento simili. Tale 


somiglianza può essere valutata mediante il calcolo di una distanza euclidea che dipende 
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dalle frequenze della struttura. Fatta questa procedura di clusterizzazione è possibile 


selezionare i poli stabili e definire un indice FAC (Frequency Assurance Criterion). 


Queste iterazioni sono ripetute successivamente fino alla definizione di un livello di soglia 


dell’ordine del modello. 


Risultati e conclusioni 


La chiesa è stata oggetto di una profonda e completa campagna di indagine. 


I risultati del monitoraggio statico, indicano uno spostamento delle torri di 0,05 m in 


direzione trasversale alla navata come mostrato in Figura 45. 


-0.06 -0.04 -0.02 0.00 
Lateral Displacement (m) 


‘0,06 -0,04 -0.02 
Lateral Displacement (m) 


Figura 45: Dettaglio delle misure degli spostamenti. 


Dalla procedura di identificazione dinamica globale è stato possibile ottenere i seguenti 


risultati riportati in Tabella 17. Le prime quattro forme modali sono state identificate con 


un errore relativo più basso del 1,5%. 


Tabella 17: Risultati dei test di vibrazione ambientale. 


Mode 1 
Mode 2 
Mode 3 
Mode 4 


CC-SSI 

f(Hz) È (3) 
2.14 1.5 
2.63 1.1 
2.85 1.4 
2.93 1.5 


FDD 


f(Hz) 


2.15 (0.5%) 
2.64 (0.4%) 
2.89 (1.4%) 
2.97 (1.4%) 
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EFDD 

f(Hz) È (3) 
2.14(<0.01%) 1.8 (20%) 
2.62 (0.4%) 1.5 (36%) 
2.89 (1.4%) 0.8 (43%) 
2.94 (0.3%) 0.8 (47%) 


A valle della procedura di model updating si è dimostrato che c'è buon accordo con i 


risultati ottenuti sperimentalmente come si può notare in Figura 46. 


Mode 1 Mode 2 Mode 3 


Experimental 
Results 
2.14 [Hz] 
FE Model Di 
= | 


2.14 [Hz] 2.55 [Hz] 2.93 [Hz] 2.94 [Hz] 


Figura 46: Confronto tra i risultati ottenuti sperimentalmente e il modello agli elementi finiti. 


Di particolare interesse sono stati i risultati del sistema di monitoraggio dinamico, che ha 
registrato, nei nove mesi in cui è stato in funzione, circa 3500 eventi. Questi sono stati 
processati automaticamente e le frequenze e i rapporti di smorzamento ricavati mostrano 
valori realistici con la variazione dei parametri ambientali. | risultati sono riportati in Figura 
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Figura 47: Evoluzione delle frequenze naturali. 
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Si può concludere quindi che il sistema di monitoraggio dinamico è stato utile per valutare 


l'evoluzione del danno e per pianificare interventi a lungo termine. 


1.2.2.6. La Chiesa del Monastero dei Girolamini a Lisbona in Portogallo 


Figura 48: Chiesa del Monastero dei Girolamini. 


Cenni storici 


La costruzione del Monastero iniziò nel 1502 e si concluse dopo circa 100 anni. In seguito 
al sisma del 1756 ci fu il collasso di una delle colonne che supportano la volta della chiesa 
e anche il parziale crollo della struttura di copertura del coro. In aggiunta a questi eventi 
durante il XIX secolo furono apportate delle modifiche alla struttura delle due torri e al 
tetto. L'effetto di questi cambiamenti nelle prestazioni sismiche della struttura rimane un 
problema aperto. Pertanto, sono state condotte negli ultimi decenni diverse indagini per 


capire meglio il comportamento strutturale. 
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Descrizione della sperimentazione effettuata 


AI fine di valutare la sicurezza della chiesa del Monastero sono state effettuate in fase 
preliminare diverse indagini in sito al fine di valutare l'integrità degli elementi strutturali, 


determinare la geometria attuale e identificare possibili vulnerabilità locali. 


È stata condotta un'analisi modale operativa per l’identificazione dinamica globale della 


struttura. In Tabella 18 si riportano le caratteristiche della strumentazione adottata. 


Tabella 18: Strumentazione adottata per l’OMA e per il Monitoraggio dinamico. 


Accelerometri triassiali 
Schede di acquisizione a 16 Bit 


Sono stati selezionati 30 punti di misura di cui 20 nella sommità dei pilastri della navata 


centrale e 10 sui muri perimetrali. 
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Figura 49: Analisi delle vibrazioni ambientali, disposizione planimetrica dei sensori. 


| dati derivanti dell'identificazione dinamica sono stati processati mediante gli algoritmi 


EFDD e SSI. 
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In questo studio, sono state condotte anche delle prove di laboratorio su dei campioni di 


materiali per definire i parametri meccanici dei principali elementi strutturali. 


Grazie alla conoscenza acquisita nella campagna di indagine preliminare, è stato realizzato 
un modello numerico ed è stata fatta l’analisi modale per la sua calibrazione. Tuttavia, a 
causa della complessità, il modello è stato aggiornato manualmente regolando il modulo 
elastico e le condizioni di vincolo. Sono state condotte diverse analisi per la valutazione 


del livello di sicurezza della chiesa tra cui l’analisi Pushover. 


Per caratterizzare completamente lo stato di salute attuale del monumento sono stati 
condotti due monitoraggi: un monitoraggio statico e uno dinamico, le cui caratteristiche 
vengono riportate in Tabella 18 e in Tabella 19. Per quanto concerne il monitoraggio 
dinamico, la strumentazione adottata è stata la medesima di quella usata per 


l’identificazione dinamica globale. 


Tabella 19: Strumentazione adottata per il monitoraggio statico. 


2 INCLINOMETRI UNIASSIALI (TL) 


6 SENSORI DI TEMPERATURA (TS) 
1DATA LOGGER 


[2 SENSORI DI MISURA DELLA TEMPERATURA DELL'ARIA E DELL'UMIDITA! 


1 ANEMOMETRO 2D 


L'analisi dei dati riguardanti il monitoraggio statico è stata effettuata per un totale di nove 
anni (da giugno 2005 a luglio 2014). L'obiettivo del monitoraggio statico era quello di 
capire se e come le condizioni ambientali influissero sulla risposta della chiesa. Per fare 


ciò sono state studiate molte correlazioni tra i parametri ambientali. 


Il monitoraggio dinamico è stato in funzione per tre anni (da agosto 2005 ad aprile 2008). 


La stima dei parametri modali viene effettuata mediante gli algoritmi EFDD e SSI. 
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Risultati e conclusioni 


In questo studio emerge chiaramente l’importanza dell'adozione di un approccio 


multidisciplinare per l’analisi strutturale e per la valutazione delle costruzioni storiche. 


Il rilievo topografico effettuato in fase preliminare ha evidenziato la presenza di 
deformazioni nelle pareti laterali di circa 6 cme il fuori piombo delle colonne che va da 5 
a 18 cm. Le indagini radar e quelle baroscopiche hanno permesso di indagare lo spessore 
delle volte che varia tra 8 e 10 cm. | test ultrasonici svolti in tre colonne e i test GPR svolti 
in altre quattro colonne della chiesa hanno rilevato una sostanziale omogeneità dei 


materiali ed un'assenza di anomalie. (Masciotta, Roque, Ramos, & Lourengo, 2016) 


Per quanto concerne i risultati ottenuti dall’identificazione dinamica globale si può notare 
un buon accordo tra quelli ottenuti con l’algoritmo EFDD e con I'SSI. In particolare, si sono 


potuti definire otto modi riportati in Tabella 20. 


Tabella 20: Risultati dell’identificazione dinamica condotta. 


Mode shape f{[Hz] E [%] MAC 
EFDD SSI EFDD SSI 
Mode 1 3.69 3.68 2.34 1.26 0.99 
Mode 2 5.12 5.04 1.11 2.68 0.92 
Mode 3 6.29 6.30 1.00 0.82 0.67 
Mode 4 7.23 7.29 0.77 1.44 0.67 
Mode 5 9.67 9.65 1.10 1.45 0.62 
Mode 6 11.64 11.65 1.20 1.46 0.36 
Mode 7 12.45 12.51 1.25 1.19 0.71 
Mode 8 14.99 15.09 1.31 Parsi 0.49 


Nonostante la buona corrispondenza tra le frequenze identificate dai due metodi, la MAC 
mostra un’ottima correlazione solo per i primi due modi. Tuttavia, questo può essere 
dovuto all'alto numero di configurazioni adottate per le misure di identificazione 
dinamica. Anche la correlazione con i fattori di smorzamento non è buona. Comunque, le 
prime due frequenze sono ben chiare, in particolare la prima forma modale è di tipo 
flessionale in direzione Nord Sud, mentre la seconda è una forma modale locale con 


significative componenti sia in direzione verticale che trasversale. 
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Quindi dopo un’intensiva campagna di indagine preliminare è stato possibile calibrare un 
modello agli elementi finiti e capire il comportamento non lineare della struttura e gli 
effetti di ciò. Di particolare interesse è quanto emerge dall'analisi pushover. Infatti, la 
risposta della chiesa è fortemente influenzata dalla colonna nord che mostra 
un'improvvisa perdita di capacità portante con conseguente aumento di spostamenti 
orizzontali e verticali fino a 13 mm e questo influenza a sua volta la risposta nella chiave 


di volta. In Figura 50, viene riportata la curva Pushover per carichi verticali. 


Load factor 


#f ——Northcolumn(top) 
Li ——South column (top) 
h Vault key 
LT = = » «North column (shaft) 
# + South column (shaft) 


-0.02 -0.015 -0.01 -0.005 0,0 0.005 0.01 0.015 0.02 


Vertical (-Z) and horizontal (Y) displacement [m] 
Figura 50: Curva di Pushover per carichi verticali. 


Sono state individuate tre frequenze corrispondenti alle forme modali 1, 3 e 4. Come 
mostrato in Figura 51 la seconda forma modale non viene riportata in quanto non è stata 


stimata accuratamente. 


Frequency Evolution over Time 
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Figura 51: Evoluzione delle frequenze durante il funzionamento del monitoraggio. 
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In Figura 52 vengono riportate le analisi condotte a valle del terremoto di magnitudo 


6.1Mw registrato nel sud ovest dell’Iberia. Come si può osservare nella Figura 52 c) i picchi 


del contenuto in frequenza dello spettro di risposta elastico rientrano esattamente nel 


range delle frequenze naturali identificate. 


Loading excitation over Time 
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Figura 52: a) RMS dei canali di acquisizione nel tempo b) risposta della chiesa durante l’evento sismico c) 
spettro di risposta elastico. 


In Tabella 21, si riporta la variazione delle frequenze tra prima e dopo il terremoto. 


Effettivamente, emerge un piccolo cambiamento, in quanto le frequenze dopo il sisma 


diminuiscono di quasi l’1% per il primo e il terzo modo. Per quanto concerne 


l’ultima 


frequenza, si nota un calo minore. Tale risultato potrebbe essere sinonimo di una piccola 


diminuzione di rigidezza della struttura. 


Tabella 21: Variazione delle frequenze prima e dopo il terremoto. 


Mode 


Ist 
3rd 
4th 


f{Hz] 


Before 


3.77 
6.40 
7.31 
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Af[®] 
After 
3.73  —-0.96 
6.34  —0.95 
7.28 —-0.53 


Si è studiata anche la correlazione frequenza-temperatura. Si osserva chiaramente che le 
frequenze crescono quasi linearmente con il crescere della temperatura e viceversa, e 


questo è un risultato attendibile e in linea con le ricerche di letteratura. 


(a) Temperture Effect on the First Mode Mode I 
* Temperature © FRQI 
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3 > > 
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Figura 53: Relazione frequenza temperatura per il primo modo. 


La variazione di temperatura ha una forte influenza sulla risposta strutturale della chiesa 


sia dal punto di vista statico che dinamico, che necessita di essere ulteriormente 


monitorata. 


LZZ:A La Cattedrale di Santa Maria Assunta a Reggio Emilia 


Figura 54: Cattedrale di Santa Maria Assunta a Reggio Emilia. 
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Cenni storici 


I primi documenti che attestano la costruzione della chiesa risalgono alla metà del IX 
secolo. Nel corso dei secoli la cattedrale è stata sottoposta ad una serie di interventi che 
hanno alterato la costruzione originale. La struttura nel corso degli anni ha subito le azioni 
di diversi sismi non ultimi quelli del 1996 (Mw = 5.4) e del 2000 (Mw = 4.5). Il quadro 
fessurativo attuale è altamente indicativo della risposta strutturale e denuncia 


chiaramente zone della chiesa che manifestano un'alta vulnerabilità sismica. 


Come si può notare in Figura 55, i danni più evidenti sono localizzati nella zona della 


cupola, nella lanterna e nella facciata che presenta anche delle fessure a taglio. 


av a CLN 


# 


Figura 55: Schema del quadro fessurativo della cattedrale. 


Descrizione della sperimentazione effettuata 


Le tecniche di indagine applicate alla Cattedrale di Reggio Emilia hanno lo scopo di 


raggiungere un alto livello di conoscenza della struttura e dei suoi elementi costitutivi. 


Le indagini disponibili da studi pregressi, come martinetti piatti singolo e doppio, 
endoscopie e campioni di provini, sono state completate attraverso l'applicazione di 


tecniche di indagine non distruttive come test sonici e test dinamici. 
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Successivamente è stata eseguita l’analisi modale operativa per definire la risposta 


dinamica della struttura. 


In Tabella 22, si riportano le caratteristiche della strumentazione adottata e in Figura 56 


si riporta la disposizione all’interno della Cattedrale. 


Tabella 22: Caratteristiche strumentazione adottata. 


Accelerometri triassiali 
Schede di acquisizione a 16 Bit 


© s.F.L ; Dc as nena 2777 norm sa ni 


c@L01 IAT 
| È i 


Figura 56: A sinistra viene proposta un’Illustrazione della posizione dei sensori nella Cattedrale di Reggio 
Emilia. | cerchi neri rappresentano i punti in cui sono stati condotti i test sonici mentre quelli in grigio 
rappresentano la posizione per i test dinamici. A destra si riporta il dettaglio della disposizione dei sensori sul 
sistema Tamburo-Cupola-Lanterna. 

La posizione dei sensori è stata preliminarmente valutata mediante un modello agli 
elementi finiti. Gli accelerometri sono collocati nei punti più rappresentativi della 
cattedrale come facciata, cupola, lanterna e la parte superiore della navata centrale. 
Particolare attenzione è stata dedicata alla cupola perché in questa zona sono presenti 


x 


diverse fessure e questo elemento è caratterizzato da un comportamento dinamico 
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complesso. Sia nella cupola che nella facciata i sensori sono stati posizionati a tre differenti 


livelli della struttura in combinazione biassiale. 


AI fine di ottenere i parametri modali della struttura i dati sono stati processati mediante 


l'algoritmo EFDD. 


È stato realizzato un modello agli elementi finiti ed è stato fatto un confronto tra le 
frequenze strutturali del modello sperimentale e i valori delle frequenze naturali 
determinate mediante l'algoritmo EFDD. Per valutare il modello numerico sulla base dei 


dati sperimentali è stata usata la MAC. 


Risultati e conclusioni 


Lo studio effettuato sottolinea l'importanza di un approccio al problema multidisciplinare. 
Infatti, il rilievo in sito è fondamentale per la definizione dei materiali e del layout 
strutturale, mentre le misure di identificazione dinamica sono altrettanto importanti per 
la determinazione dei parametri modali della struttura e quindi per la messa a punto di 


un modello agli elementi finiti. (Casarin & Modena, 2008) 


Tabella 23: Test di identificazione dinamica, frequenze sperimentali EXP, frequenze FEM, Massa partecipante 
e Descrizione modale. 


Frequency Frequency Mass part. (%) 

Mode EXP Detection FEM 

nr. (Hz) area (setups) (Hz) Longitudinal Transversal Mode description 

1 2.64 Dome 2.65 0.21 16.24 Transversal bending, dome 

2 2.68 Dome 2.66 7.76 0.28 Longitudinal bending, dome 

3 2.97 Dome 2.90 1.36 20.05 Transverse bending, dome 

2 3.05 Facade 3.02 27.82 1.20 Longitudinal bending, facade 

b) 3.12-3.14 Facade/Nave 3.33 0.02 9.05 Transverse bending, nave 

6 3.79 Facade 3.89 16.17 1.07 Apses area: longitudinal 

bending, facade 

LI 4.08 3.67 0.01 Global bending 

8 4.04 Nave 4.36 4.11 12.75 Transverse bending dome 
torsion. facade 

9 4.59 Nave 4.51 2.00 4.64 Transverse bending dome 
torsion. facade 

10 4.63 0.36 2.03 Torsion, facade 


Dai risultati riportati in Tabella 23, si nota che solo le prime quattro frequenza modali 


hanno un buon accordo e valori di MAC prossimi all'unità. 
p 
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2.IL CASO STUDIO: LA CHIESA DI SANTA MARIA 
IN VIA IN CAMERINO (MC) 


In questo capitolo, si descrive il caso studio del presente lavoro di Tesi. 


Inizialmente viene fatta una panoramica sulla storia della chiesa di Santa Maria in Via, 
andando a descrivere nel dettaglio sia gli eventi che l’hanno interessata nel corso della 
storia sia gli aspetti architettonici. Viene poi descritto lo stato di danno relativo al sisma 
del Centro Italia del 2016 che ha portato poi nel 2017 all’installazione di un sistema di 


messa in sicurezza provvisoria. 


2.1 Inquadramento storico e conformazione architettonica 


La città di Camerino si trova in una suggestiva posizione al centro della zona montana della 
provincia di Macerata è situata tra le valli del Chienti e del Potenza è chiusa a sud dal 


massiccio dei Monti Sibillini e a nord dal Monte San Vicino. 


In questo paesaggio, si trova nel centro storico di Camerino, la Chiesa di Santa Maria in 


Via che si affaccia sulla piazza Umberto |. 


Le prime notizie dell’esistenza dell’edificio di Santa Maria in Via sono riprodotte nelle 
“Rationes Decimanorum” del 1299 e in un testamento del 1341, nove anni prima del 
ritorno dei mille crociati che secondo la tradizione riportarono da Smirne la venerata 
immagine che raffigura la Madonna col bambino. Nel XVI secolo la chiesa era una 
costruzione povera ed irregolare. L'attuale edificio è stato costruito tra il 1639 e il 1642 
sull'area precedente ai cui lavori cui partecipò Camillo Arcucci da Sigillo, figura molto 
vicina a Francesco Borromini e venne consacrata nel 1654 dall'allora vescovo di Camerino 


Monsignor Altieri. (Petrucci) 
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La pianta è caratterizzata da andamento trapezoidale ed è frutto dell’accorpamento di 
alcuni edifici preesistenti che vennero uniti al fine di ampliare e regolarizzare la chiesa 


originale. 


Il terremoto del 28 luglio 1799 provocò ingenti danni alla chiesa. La copertura crollata fu 
sostituita da un tetto a doppia falda sostenuto da capriate lignee e completata da una 
finta cupola in camorcanna. La cupola al colmo presenta un ovale che interrompe e 


raccorda le lesene che ne suddividono la superficie interna in otto spicchi. 


Figura 57: La cupola in camorcanna prima del crollo. 


Attualmente la chiesa è costituita planimetricamente da un unico corpo di fabbrica dal 
quale si elevano il tiburio, a pianta ottagonale irregolare, il campanile e il corpo della 


facciata. 


L’interno della chiesa presenta sostanzialmente l’aspetto architettonico primitivo, 
raccolto ed armonioso e le decorazioni del secolo scorso. La pianta ha forma ellittica con 
presbiterio absidato sull'asse maggiore e quattro cappelle radiali a nicchia emiciclica 
ricavate nello spessore della muratura. Sopra agli accessi laterali, in coincidenza dell’asse 


minore dell’ellisse, sono posti due coretti. La struttura perimetrale interna è scandita da 
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paraste che sorreggono una trabeazione a fregio continuo su cui si imposta la volta in 
camorcanna. Il tiburio che si eleva per oltre 8 m rispetto alla cornice interna è 
caratterizzato dalla presenza di quattro finestroni in corrispondenza degli altari laterali, 
alcune testimonianze suggeriscono che in antichità i finestroni sarebbero stati sei due dei 
quali tamponati al giorno d’oggi. All'esterno è dotato di contrafforti in corrispondenza 


degli angoli. (Petrucci) 


II corpo della facciata può essere diviso in due ordini principali: l'ordine inferiore è 
inglobato nella pianta trapezoidale della chiesa, mentre l'ordine superiore si eleva in 
adiacenza al tiburio per un’altezza di 23 m. Inoltre, il corpo della facciata ha pianta 
trapezoidale isoscele con lunghezza media delle basi di circa 16 m in direzione trasversale 
della chiesa ed altezza di circa 6 m in direzione longitudinale rispetto alla pianta 
dell’edificio. L’interno del corpo di facciata è caratterizzato da tre orizzontamenti in legno 
in cui il primo si trova a circa 4,5 m dal piano di calpestio dell’aula, il secondo a 10,35 m e 
il terzo alla quota dell'orologio a circa 15,5 m. In pianta il corpo di facciata è costituito da 
tre vani a ciascuna delle elevazioni separati da due pareti di spina. Nel blocco superiore si 
nota la presenza di un sistema di catene in direzione longitudinale ancorate con capi 
chiave a paletto all’esterno. All’interno della parete di facciata, sul lato sinistro è presente 
una scala a chiocciola che consente l’accesso alla cantoria e al piano superiore. La 
presenza di questo elemento produce una significativa discontinuità nella muratura. 

In merito al corpo murario il paramento esterno è interamente a faccia a vista e risulta 


essere realizzato con una muratura molto regolare e compatta di mattoni e malta di calce. 


Il corpo della facciata invece è costituito da pietra rivestita con laterizio di spessore molto 


diverso. 
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Figura 58: Facciata della Chiesa di Santa Maria in Via. 


Il campanile della chiesa ha pianta ottagonale e dall’osservazione delle caratteristiche 


architettoniche, la sua costruzione potrebbe coincidere con quella della chiesa. 


La Chiesa ha subito danni diffusi in seguito all'evento del 1997 (terremoto Umbria - 


Marche) e dopo un lungo restauro venne riaperta al culto nel 2006. 


Solo dieci anni più tardi, a causa del terremoto del Centro Italia del 2016, l’edificio è 


nuovamente inagibile. 
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2.2 Lo stato di danno 


2.2.1 Il terremoto del Centro Italia del 2016 


Nel 2016 una serie di eventi sismici ha interessato il Centro Italia provocando morti e 


distruzione. 


La prima forte scossa è avvenuta il 24 agosto 2016 alle ore 01:36:32 GMT, e ha avuto una 
magnitudo di 6.0 con epicentro situato lungo la valle del Tronto, tra i comuni di Arquata 


del Tronto (AP), Accumuli (RI) e Amatrice (RI) causando 300 morti. 


Il 26 ottobre 2016 sono avvenute due potenti repliche di magnitudo 5.4 (17:10:36 UTC) e 
5.9 (19:18:06 UTC) che hanno esteso il volume sismogenetico verso Nord-Ovest. Quattro 
giorni più tardi il 30 ottobre 2016 un nuovo evento con magnitudo 6.5 ha colpito la zona 


con epicentro a Norcia (PG). (lervolino, Baltzopoulos, Chioccarelli, & Suzuki, 2016) 


Il 18 gennaio 2017 si sono verificati altri quattro eventi con magnitudo superiore a 5 ed 


epicentri localizzati tra i comuni aquilani di Montereale, Capitignano e Cagnano Amiterno. 


Nel periodo agosto 2016 — agosto 2020 la rete sismica nazionale ha registrato oltre 118 
mila eventi sismici di cui 1182 hanno avuto magnitudo compresa tra 3 e 3.9, 66 tra 4 e 
4.9, mentre 9 eventi hanno avuto magnitudo compresa tra 5 e 5.9 e 2 hanno superato 
magnitudo 6. In Figura 59 si riporta la mappa degli eventi registrati nella zona nei quattro 


anni di osservazione. 


La sequenza sismica ha colpito un’area dove in passato si sono verificati diversi grandi 
terremoti. Infatti, è ben documentata una sequenza sismica che ha colpito l'Appennino 
centro meridionale nel XVIII secolo. Il terremoto con maggiore intensità che interessò la 


zona si verificò nel 1703 e distrusse diversi villaggi nelle vicinanze di Norcia (PG). 
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@ INGV terremoti 


SEQUENZA SISMICA IN 
ITALIA CENTRALE 


Sismicità dal 24 agosto 2019 
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® da20a29 
@ 0da30239 
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fonte dati: http/terremoti.ingv.it 


Figura 59: Sequenza sismica agosto 2016 — agosto 2020. 


2.2.2 Lo stato di conservazione e di danno 


Prima di subire i danni relativi al sisma del 2016 la costruzione si presentava in ottimo 
stato di conservazione in quanto restaurata dopo l'evento del 1997. Tali interventi 
coinvolsero maggiormente la volta ellittica in camorcanna che a causa della sua leggerezza 
riportò le maggiori conseguenze. Vennero effettuati interventi di consolidamento della 
volta e vennero rafforzati gli ancoraggi originali tra i travetti e le centine lignee. (Marche, 


2007) 


I primi danni nella chiesa si sono registrati dopo la scossa del 24 agosto 2016 che ha 


portato a una chiusura dell’edificio. 
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In seguito alla seconda scossa del 26 ottobre 2016 anche il campanile collassò (Figura 60), 
mentre il basamento superstite che si eleva per circa quattro metri, non fu interessato in 


quanto sono presenti due ordini di catene che hanno evitato il collasso totale. 


Figura 60: Resti del campanile dopo le scosse del 26 ottobre 2016. 


In conseguenza degli eventi verificatesi nel gennaio 2017, in aggiunta alle pessime 
condizioni atmosferiche, l’edificio ha esibito ulteriori danni che hanno interessato 
principalmente la cupola con il collasso di una parte di essa, il crollo parziale della 
copertura lignea nonché lesioni diffuse sulla facciata e sulle pareti interne ed esterne 


come si può vedere in Figura 61. 


68 


Figura 61: Veduta aerea dei crolli in copertura. 


Come riportato nella relazione tecnica illustrativa, (Leoni & Gioiella, 2017), notevoli danni 
si evidenziano all'altezza del tiburio, in particolare nella porzione posteriore che è crollata 
verso l’interno dell’aula in seguito ai progressivi eventi sismici e per le forti nevicate del 


gennaio 2017. 


Le strutture verticali del tiburio sono costituite prevalentemente da muratura in pietra di 
blocchi sbozzati caratterizzata da inserti di blocchi in pietra non lavorati e laterizio. Il crollo 
della cupola è in parte dovuto alla pessima qualità dei materiali costituenti ed in parte alla 
presenza di tre puntoni diagonali sul padiglione della copertura. In questa zona si nota 
anche la presenza di alcune nicchie tamponate a dimostrazione del fatto che la chiesa ha 
subito modifiche nel corso dei secoli, tuttavia ciò ha aumentato la presenza di punti di 
discontinuità nella muratura producendo un ulteriore indebolimento della costruzione. In 
questa zona in seguito al sisma del 1997 è stato eseguito un intervento di cerchiaggio 
costituito da un piatto d'acciaio connesso con barre verticali immerse nella muratura, ma 
ciò non è stato sufficiente ad evitare il crollo. Nel 2017 è stato eseguito prima un intonaco 


nella zona della muratura crollata per evitare che il dilavamento delle superfici possa 
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incrementarne la vulnerabilità successivamente è stata realizzata una struttura 


provvisoria al fine di proteggere l'interno dell'aula. 


Figura 62: Dettagli dei meccanismi di danno. 


Le pareti interne ed esterne al corpo della navata centrale e la parte absidale non sono 
interessate da quadri fessurativi significativi, ad eccezione dei due altari laterali in 
adiacenza del corpo di facciata che hanno subito importanti danni agli archi e alla volta. Il 
cornicione interno è stato danneggiato dalla caduta di blocchi dall'alto e presenta alcune 


porzioni pericolanti. (Leoni & Gioiella, 2017) 


Le pareti laterali del corpo dell’edificio e della parte absidale non sono interessate da crolli 
o da quadri fessurativi importanti eccezion fatta per la cella campanaria dove sono 
presenti alcune lesioni diagonali che sono focalizzate nelle fasce di piano tra un’apertura 


e l’altra. 
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PROSPETTO EST PROSPETTO NORD 


Figura 63: Quadro fessurativo prospetti Est - Nord — Ovest — Sud. 


La facciata presenta un meccanismo di ribaltamento fuori dal piano che è parzialmente 
pregresso ma ha subito un aggravio in seguito agli eventi sismici del 2016, andando ad 


interessare entrambi gli ordini superiore e inferiore. In seguito al meccanismo descritto 
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questo elemento ha subito una rotazione esterna al piano nell’ordine inferiore, mentre la 
parte superiore è stata oggetto di una traslazione rigida la quale ha interessato anche una 
parte del tiburio adiacente. Tutto ciò ha portato anche alla formazione di diverse fessure 
nel rivestimento in laterizio generando anche diverse disgregazioni. Attraverso un’attenta 
analisi si è potuto osservare anche il distacco di una parte del cornicione al primo ordine 
il quale sembra non essere ammorsato alla parete. Anche il portale ha subito dislocazioni 
significative con alcuni cedimenti dell’architrave. Vale la pena sottolineare anche che il 


corpo della facciata presenta un fuori piombo di 30 cm. (Leoni & Gioiella, 2017) 


II meccanismo di ribaltamento fuori dal piano accompagnato dal meccanismo di 
traslazione rigida della parte superiore ha provocato danni significativi soprattutto nella 
parte destra della facciata e generando anche diversi meccanismi secondari di 


ribaltamento fuori dal piano di porzioni di pareti laterali. 


Figura 64: Dettagli del meccanismo di danno presente sulla facciata. 
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2.3 Il presidio di messa in sicurezza 


Nei mesi successivi agli eventi che hanno colpito il territorio, alla luce dello stato di danno 
descritto nella relazione tecnica illustrativa, (Leoni & Gioiella, 2017), è stato messo a 
punto un presidio provvisorio con lo scopo di mettere in sicurezza le vie d'accesso alla 
città e le aree pubbliche adiacenti alla chiesa e ridurre al minimo l'occupazione di spazi 
per il traffico. Sono stati realizzati inoltre, sia interventi di bonifica e stabilizzazione delle 


murature esterne che di protezione dello spazio interno della chiesa. 


L’opera che è stata realizzata ha come scopo quello di evitare ulteriormente l’attivazione 
dei meccanismi del corpo di facciata, del tiburio nonché eventuali meccanismi secondari 


non prevedibili ed eventuali crolli che potrebbero scaturirne. 


Gli interventi principali hanno previsto l'adozione di un sistema di incatenamento e 
trattenimento del corpo di facciata e della porzione di tiburio adiacente. Inoltre, alla luce 
dell'importante fuori piombo, il sistema adottato ha anche lo scopo di non far progredire 


questa criticità. 


Per il caso studio in esame viste le dimensioni geometriche e le caratteristiche 
morfologiche è stata adottata una soluzione non convenzionale per la messa in sicurezza. 
È stato inserito un sistema di tubi-giunti all'interno dell'aula e nella parte sommitale 
esterna del tiburio, che ha consentito la stabilizzazione e la successiva ricostruzione della 
porzione di muratura crollata del tiburio, nonché la realizzazione di una copertura 
provvisoria, al fine di impedire ulteriori infiltrazioni di acque meteoriche all'interno 


dell'aula. 


Il presidio che è stato realizzato nella facciata come mostrato in Figura 65, è costituito da 
un sistema a graticcio, composto da tre ordini di elementi strutturali le cui estremità sono 
collegate a due controventi in acciaio disposti lateralmente. Questo sistema di tralicci e 
controventi è stato fondato superficialmente attraverso travi in calcestruzzo armato. Esso 
è trattenuto ai lati da due ordini di tiranti longitudinali costituiti da sei funi a trefolo di 


diametro 26 mm e da otto funi a trefolo con diametro 26 mm, per ciascuno dei due lati, 
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sono stati predisposti due ordini di tiranti trasversali costituiti da due funi a trefolo di 
diametro 26 mm ciascuno. Il graticcio è costituito da quattro travi reticolari 
tridimensionali in acciaio che si trovano sul fronte principale della chiesa e sono collegate 
ai due controventi reticolari laterali e solo le prime due travi sono collegati ai tiranti 


longitudinali. (Leoni & Gioiella, 2017) 
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Figura 65: Elementi di presidio prospetto frontale della facciata. 


| tiranti longitudinali sono stati messi esternamente lungo il perimetro della chiesa con un 
contrasto sul prospetto posteriore offerto da un sistema di deviatori. Viceversa, i tiranti 
trasversali attraversano la struttura allo scopo di collegare mutuamente le ultime due 
elevazioni dei controventi reticolari per cercare di evitare il ribaltamento delle parti 


laterali della facciata. 
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Figura 66: Elementi di presidio prospetto sinistro. 
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Figura 67: Elementi di presidio prospetto posteriore. 


75 


Gli elementi secondari sono costituti da due livelli di profili verticali in acciaio, disposti sia 
in facciata, che lateralmente al di sopra e al di sotto del cornicione. Il terzo ordine degli 
elementi che compongono il graticcio è rappresentato da pannelli in legno tipo XLAM 
posti tra i profili verticali e la facciata (principale e laterale) della chiesa utilizzati in quanto 
più adattabili alle superfici in muratura danneggiate e nel caso di ulteriori scosse il loro 


utilizzo impedisce la caduta di ulteriore materiale. 


Le azioni orizzontali che si sviluppano sulle pareti laterali e sulla parete posteriore, in 
conseguenza all'azione di tiro dei cavi longitudinali, possono trovare dei punti di debolezza 


nelle nicchie e nelle aperture delle pareti. 
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Figura 68: Presidio di sicurezza adottato. 
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3. IL SISTEMA DI MONITORAGGIO DELLA 
CHIESA DI SANTA MARIA IN VIA 


AI fine di prevenire un'ulteriore evoluzione dello stato di danno descritto, oltre 
all’introduzione dell'imponente sistema di messa in sicurezza in parallelo è stato messo a 
punto un piano di monitoraggio per il controllo dell'evoluzione nel tempo di alcune 


proprietà statiche e dinamiche della chiesa e del presidio. 


Nell'ottobre del 2020 è stata eseguita una prima campagna di indagine sperimentale 
consistente in prove di vibrazione ambientale su tutta la chiesa e dal novembre del 2020 
fino al 6 gennaio del 2021 sono stati installati i primi sensori del sistema di monitoraggio 


permanente per la valutazione dello stato di salute della facciata. 


Nel caso studio proposto le procedure di identificazione dinamica e il monitoraggio dello 
stato di salute SHM vengono condotte attraverso la registrazione di storie temporali di 
accelerazioni mediante sensori elettrici. Infatti, si acquisiscono le risposte dinamiche della 


struttura sottoposta a sorgenti di vibrazione naturale o di origine antropica. 


Una volta effettuata l'acquisizione dei segnali, al fine di ottenere le caratteristiche modali 
(frequenze, forme modali, smorzamenti), viene applicato un algoritmo di identificazione 
automatico che verrà descritto nel seguito. Questo è stato possibile attraverso l'utilizzo 


del software di calcolo Matlab. 


Nel presente capitolo vengono quindi mostrati i principi teorici alla base delle due 
campagne di indagine svolte in accordo con le indicazioni fornite da Farrar e Worden i 
quali individuano quattro step per la definizione di un sistema di Structural Health 


Monitoring: (Farrar & Worden, 2013) 
e Valutazione operativa; 
e Acquisizione dei dati; 


e Selezione delle caratteristiche; 
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e Applicazione dei modelli statistici per la validazione dei risultati. 


Negli step 2 e 4 sono incluse operazioni di normalizzazione, pulizia compressione e fusione 
dei dati. Questi approcci sono stati implementati attraverso una combinazione di sistemi 


hardware e software. 


SL, Identificazione dinamica preliminare 


Nella Chiesa di Santa Maria in Via è stata condotta una campagna di indagine dinamica 


preliminare all’installazione del monitoraggio dinamico dello stato di salute della chiesa. 


L’identificazione della dinamica globale della chiesa è stata condotta attraverso prove di 
vibrazione ambientale attraverso l’uso di 16 accelerometri piezoelettrici disposti secondo 


due configurazioni mostrate in Figura 69. 


CONFIGURAZIONE 1 CONFIGURAZIONE 2 


Figura 69: Disposizione degli acceleromettri. 


Per ottenere informazioni sul comportamento globale della chiesa sono stati monitorati 


otto punti situati intorno all’ellisse del tiburio a due diverse altezze dal livello del suolo. 
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In ogni punto (da A a H) sono stati posizionati due accelerometri monoassiali con direzioni 


di misura lungo gli assi principali della sala della chiesa (chiamate direzioni x e y). 
Nel dettaglio: 


e Configurazione 1: gli accelerometri sono stati posizionati in otto punti diversi (da 


1A a 1H) a 12m daterra. 


e Configurazione 2: due accelerometri sono stati lasciati nella stessa posizione (1B) 
per essere usati come riferimento mentre gli altri sono stati spostati a 19,7 m da 


terra; le posizioni 2C e 2H sono state escluse in quanto inaccessibili. 
La strumentazione adottata è consistita in: 


e Sedici accelerometri PCB 393B31 con una sensibilità di (+5%) 10V/g; Broadband 
Resolution: 0.000001 g rms; Range di misura: 0.5 g pk; Range di frequenza: (+5%) 
da 0.1 a 200 Hz: 


Figura 70: Accelerometro PCB 393B31. 


e Cinque schede di acquisizione NI9234 a quattro canali per la conversione del 
segnale A/D, con frequenza di risoluzione 24 bit, range del segnale +5V e 


frequenza di campionamento massima di 51,2 kS/s/ch.: 


80 


bs iena 
SV2680.00 
EITtete alitizanai 
ita 


Figura 71: Scheda di acquisizione NI9234. 


Una CompactRIO 9045, con 8 slot; processore Dual Core 1.30 GHz; 2 GB DRAM; 4 
GB Memoria; e temperatura di esercizio da -20°C a 55°C, Kintex-7 70T FPGA. 


Figura 72: CompactRIO 9045. 


Tre CompactDAQ 9185 a quattro slot. 


Figura 73: CompactDAQ 9185. 
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È stata creata una rete di sensori di misura. Nel dettaglio è stata posizionata la 
CompactRIO nei pressi della stazione di misura di riferimento 1B, e tre CompactDAQ vicino 
ai restanti punti di misura. La sincronizzazione è stata ottenuta mediante la tecnologia 
TSN che fornisce una sincronizzazione temporale distribuita e una comunicazione 


deterministica utilizzando reti Ethernet standard. 


Attraverso le schede di acquisizione NI9234 il segnale è stato campionato a 2048 Hz. Prima 
dell’identificazione dei parametri modali il segnale è stato pre-processato; in particolare 
sono stati eliminati i trend spuri sottraendo al segnale un polinomio di terzo grado; in 
seguito il segnale è stato filtrato con un filtro passa basso con una frequenza di taglio di 
25 Hz e ricampionato a 51,2 Hz al fine di ridurre la quantità di dati e rendere l’analisi più 


veloce. 


L’identificazione dei parametri modali è stata effettuata utilizzando la tecnica SSI-PC 


descritta al paragrafo 3.4.1.3. 


Poiché le acquisizioni relative alle due configurazioni sono asincrone, al fine di ottenere 
una corretta forma modale, è stato necessario procedere con un'operazione di marge dei 
dati; in questo caso sono riportati i risultati ottenuti con la tecnica PoSER (Post Separate 


Estimation Re-scaling). 


In Figura 74 vengono mostrati i diagrammi di stabilizzazione ottenuti dalle prove in 


Configurazione 1 e Configurazione 2. 

| criteri adottati per definire un modo “stabile” sono stati: 
e Af<z1% 
e MAC> 95% 


Nella definizione del diagramma di stabilizzazione è stato riportato solo il primo valore 


singolare SV1 e i modi stabili sono indicati con un cerchio di colore blu scuro. 
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Nella stessa figura, sono mostrati anche i diagrammi frequenza-smorzamento, dove si 


possono apprezzare i valori dei rapporti di smorzamento relativi alle soluzioni stabili. 


CONFIGURAZIONE 1 CONFIGURAZIONE 2 
100 
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Figura 74: Diagrammi di stabilizzazione per le due configurazioni studiate. 


Da questa identificazione dinamica modale sono stati individuati principalmente sei modi 


di vibrare le cui caratteristiche sono riportate in Figura 75. 


Vale la pena notare che nonostante il quadro fessurativo cui è soggetta la chiesa, si 
riescono ad individuare i primi due modi traslazionali che sono disaccoppiati e il terzo 


modo torsionale. Questo è evidente anche dai valori dalla matrice Auto-Mac riportata in 


Figura 76. 
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1° Modo 2° Modo 3° Modo 
f= 2,13 Hz f= 2,62 Hz f= 3,33 Hz 
£= 0,81% €= 0,86% €£= 0,85% 


4° Modo 5° Modo 6° Modo 
f = 3,66 Hz f= 3,96 Hz f= 4,36 Hz 
€= 1,11% £= 0,69% €£= 1,01% 


ivi 


Figura 75: Forme modali individuate dall’identificazione dinamica globale. 


ceo] os [n] 00 DB 


F[Hz] 213 262 333 366 396 436 


Figura 76: Matrice Auto-Mac per l’identificazione dinamica globale della struttura. 
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Per quanto riguarda i modi superiori al terzo, questi sembrano coinvolgere la dinamica 
locale del tiburio e della facciata. Infine, dalla Figura 74 si può notare che sono presenti 
modi superiori che, in questa fase, non sono stati riportati in quanto necessitano 


un'ulteriore investigazione 


GIA Valutazione Operativa 


Il processo di valutazione operativa rappresenta il primo step da compiere nel caso si 
voglia introdurre un sistema di Structural Health Monitoring secondo le indicazioni di 


Farrar & Worden. (Farrar & Worden, 2013) 


In seguito all’identificazione dinamica eseguita, sono stati individuati i primi sei modi di 
vibrare dell’edificio, tuttavia diverse sono le ragioni che hanno portato all’installazione del 


sistema SHM. 


Lo scopo principale del sistema di monitoraggio che è stato installato nel novembre del 
2020 è quello di monitorare l'evoluzione del meccanismo di ribaltamento della facciata 
attivato in seguito agli eventi sismici del 2016 e che rappresenta una delle maggiori 
criticità della chiesa. Conseguentemente a ciò si vuole cercare di valutare e monitorare 


l'eventuale interazione del sistema di messa in sicurezza con la struttura. 


Il secondo motivo non meno importante che giustifica l'introduzione di un sistema di 
monitoraggio risiede nel fatto che la Chiesa di Santa Maria in Via rappresenta un 
patrimonio artistico e culturale di notevole importanza nelle Marche e in Italia che deve 
essere conservato. Pertanto, anche alla luce del grave stato di danno individuato, è 
importante valutare l’evoluzione dei parametri modali nel tempo ed eventuali 


correlazioni che essi possono manifestare con i parametri ambientali. 


CICA Acquisizione Normalizzazione e Pulizia dei dati 


Per effettuare queste analisi di acquisizione dei dati si è adottata la metodologia OMA 


“Operational Modal Analysis” o "Analisi modale operativa" che consente di determinare 
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le caratteristiche modali di una struttura, intese come frequenze smorzamenti e forme 
modali, in condizioni operative. Questa tecnica si basa su un approccio di tipo “output- 
only” in quanto si misura sulla struttura solo la risposta senza che vengano applicate 
forzanti esterne. L’eccitazione è di natura casuale e dipende dal contesto in cui essa si 
trova. Le sorgenti di vibrazione possono essere sia naturali che artificiali (vento, sisma, 


traffico, rumore di macchine operative ecc...). 


La fonte di eccitazione non viene misurata, ma si conoscono a priori alcune caratteristiche 
spettrali di cui la più importante è quella di avere uno spettro in frequenza costante che 
viene chiamato rumore bianco e quindi analizzando la risposta del sistema si dovrebbe 
pervenire ai parametri modali. In Figura 77 si riporta un diagramma relativo al 


funzionamento di un’analisi di tipo OMA. 


Due ipotesi vengono fatte alla base dell’analisi modale operativa ovvero che il sistema si 
comporta in modo lineare e invariante nel tempo. La linearità implica che se l'input ha 
una determinata ampiezza genera un output con la medesima ampiezza. L’invarianza nel 


tempo implica che il sistema fisico non cambi nel tempo. 


Questa tecnica è molto utile in quanto consente di ottenere: 


Caratterizzazione dinamica della struttura attraverso determinati algoritmi nel 


dominio del tempo o nel dominio delle frequenze; 
e Velocità nell'esecuzione delle prove; 
e Economicità delle prove; 


e  Molteplicità dei campi di applicazione (Identificazione parametri modali, 


Monitoraggio strutturale e Collaudo dinamico di nuove strutture). 
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SISTEMA DINAMICO 


OUTPUT 
INPUT Caratteristiche dinamiche 
Rumore bianco fi- d- È 


Figura 77: Diagramma di flusso analisi modale operativa. 


3.3.1. Sistema di acquisizione nel sistema SHM 


Il sistema di acquisizione che è stato adottato per il monitoraggio strutturale in opera dal 


16 novembre 2020 al 6 gennaio 2021 è consistito in: 


e Dueaccelerometri piezoelettrici di tipo PCB 393A03 disposti ortogonalmente alla 
facciata all'altezza di 12,6 m da terra e collegati a schede di acquisizione di tipo 
N19234 le cui caratteristiche principali sono: Sensibilità (+5%) 1.0 V/g; 
Broadband Resolution: 0.00001 g rms; Range di frequenza (+5%) da 0.5 a 2000 
Hz; Range di misura 5 g pk. In Figura 78 nell'immagine b si mostra il sensore sulla 


facciata. 


a) b) 


Figura 78: a) Dettaglio accelerometro PCB 393403; b) applicazione alla facciata della Chiesa di Santa Maria 
in Via 
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e Unascheda di acquisizione NI 9234 a quattro canali per la conversione del segnale 
A/D, con frequenza di risoluzione 24 bit, range del segnale +5V e frequenza di 


campionamento massima di 51,2 kS/s/ch.: 


Ao) 
< 
Nice 
SERIA i Ci 
e e 


Figura 79: Scheda di acquisizione NI9234 


e Due Termocoppie di T con terminale a piastrina di cui una disposta internamente 
alla struttura e una esternamente sul presidio in acciaio. Sono collegate a una 
scheda di acquisizione di tipo NI9210. In Figura 80, si riporta un dettaglio della 


termocoppia all’interno dell’aula dell’edifico. 


a) 


Figura 80: a) Termocoppia di tipo T con terminale a piastrina b) dettaglio della termocoppia all’interno 
dell’aula della Chiesa. 
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e Una scheda di acquisizione NI9210 a quattro canali con doppio modulo di 
acquisizione, con frequenza di risoluzione a 24 bits, range del segnale + 80mV C e 


frequenza di campionamento 14 S/s. 


N 


» 
|» 
} 


sì 
» 


Figura 81: Scheda di acquisizione 9210 


e Controller CompactDAQ-9133 a otto slot, dual-core Atom 1.33 GHz, con 16GB di 


memoria. 


Figura 82: CompactDAQ-9133 


e Gruppo di continuità UPS per garantire l'acquisizione dei dati anche in presenza 


di temporanei guasti elettrici. 
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Figura 83: Gruppo di continuità UPS 


e Connessione Internet. 
Una rappresentazione grafica del sistema di monitoraggio è riportata in Figura 84. 


| sensori sono stati disposti secondo questo schema perché lo scopo primario del sistema 
di monitoraggio è quello di monitorare il meccanismo di ribaltamento della facciata. Come 
conseguenza di questa disposizione degli accelerometri, il primo modo individuato 
nell’identificazione dinamica non potrà essere rilevato in quanto esso non coinvolge la 
facciata. Si specifica che si tratta di una configurazione preliminare: i risultati 
dell’identificazione preliminare, uniti ai primi risultati del monitoraggio, consentiranno di 


definire un miglior posizionamento dei sensori e un ottimale numero dei punti di misura. 
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Figura 84: Schema di disposizione del sistema di monitoraggio 


II software di acquisizione implementato nella CDAQ è stato sviluppato in ambiente 
Labview e acquisisce dati in continuo. 

Sono stati scelti accelerometri non troppo sensibili perchè si è voluta eseguire 
un’acquisizione a tutto tondo cercando di registrare sia le vibrazioni ambientali sia 
eventuali sismi. AI contrario se si fossero usati accelerometri con maggiore sensibilità in 
caso di evento sismico si sarebbe raggiunta facilmente la saturazione perdendo quindi la 
registrazione della risposta strutturale. 

Il sistema cosi composto ha consentito la trasmissione e l'archiviazione dei dati acquisiti 
dal sistema SHM mediante protocollo FTP (File Transfer Protocol) ad un server che si trova 
presso il dipartimento DICEA dell’Università Politecnica delle Marche secondo lo schema 


riportato in Figura 85. 
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ACQUISIZIONE DATI \ É SERVER NY 
Chiesa di Santa Maria in Via - | FTP Università Politecnica delle Marche 
Camerino DICEA - Ancona 
LL ) L Ian] È 


Figura 85: Schema trasmissione Dati 


Questo approccio di trasmissione dei dati in continuo permette di indagare in tempi brevi 
eventuali dipendenze dalle condizioni ambientali, ma permette anche di registrare 


eventuali eventi sismici senza il rischio di perdere dei campioni di registrazioni. 


3.3.2. Normalizzazione dei dati registrati 


Il processo di normalizzazione consiste nell’interpretazione dei dati derivanti dal sistema 
di acquisizione e in particolare consiste nel separare le modifiche introdotte all’interno 
della registrazione causate da un eventuale stato di danno rispetto a quelle causate dalle 
variabili ambientali. 

Come accennato in precedenza, la struttura oggetto di studio essendo molto danneggiata 
ha un comportamento fortemente non lineare, inoltre la chiesa trovandosi all’interno 
della zona rossa non è soggetta a vibrazioni ambientali di origine antropica cosicché 
eventuali variazioni del segnale misurato devono essere valutate attentamente. 

Alla luce di ciò, la normalizzazione dei dati è molto importante e senza questa operazione 
si rischierebbe di confondere gli effetti provocati dalle vibrazioni ambientali dagli effetti 
prodotti da un nuovo stato di danno. Questa fase, verrà applicata a valle dei periodi di 
addestramento, quando saranno state indagate le dipendenze dei parametri monitorati 


dalle condizioni ambientali. 


3.3.3. Pulizia dei dati registrati 


La pulizia dei dati è un processo di scelta selettiva dei dati da cui poi estrarre le 


caratteristiche modali. 
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Mediante le schede di acquisizione NI9234 il segnale è stato campionato alla frequenza di 
2048 Hz. È stato ricampionato alla frequenza di 51,2Hz ed è stato applicato un filtro 
passabasso alla frequenza di Nyquist di 25,6 Hz al fine di ridurre la quantità di dati e 


rendere l’analisi più veloce. 


Le time history come già accennato al paragrafo 3.3.1 sono state acquisite in continuo in 
finestre temporali di dieci minuti. AI fine di effettuare l’identificazione dinamica su una 
time history sufficientemente lunga, le finestre temporali sono state accoppiate e 


l'estrazione delle caratteristiche modali è stata condotta in time history di venti minuti. 


3.4. Estrazione delle caratteristiche 


Una delle operazioni più importanti nell’analisi modale operativa è l’implementazione di 


algoritmi finalizzati all'estrazione delle caratteristiche modali della struttura cioè: 


e Frequenze; 
e Formemodali; 


e Rapporti di smorzamento. 


Questi parametri rappresentano una carta d'identità della struttura perché essi 


dipendono da altre caratteristiche come massa e rigidezza. 


Gli algoritmi che possono essere adottati per l’identificazione di tali parametri si basano 


su due tecniche principali: 


e Tecniche nel dominio del tempo come ad esempio: Decremento logaritmo; LSCE 
(Least Squares Complex Exponential); SSI- DD (Stochastic Subspace Identification 
Data Driven); SSI — COV (Stochastic Subspace Identification Covariance); SSI — PC 


(Stochastic Subspace Identification Principal Component); 


e Tecniche nel dominio delle frequenze come ad esempio il Peak Picking (BFD Basic 
Frequency Domain); LSCF (Least Squares Complex Frequency); FDD (Frequency 


Domain Decomposition); EFDD (Enhanced Frequency Domain Decomposition). 
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Per il caso studio proposto, si è scelto di utilizzare l'algoritmo SSI-PC, in quanto le tecniche 
di estrazione dei parametri modali basate sul dominio del tempo sono più sensibili e 


riescono a cogliere bene anche frequenze molto vicine tra loro. 


Densità di potenza spettrale e decomposizione dei valori singolari 


A supporto dell’identificazione dinamica vengono usualmente utilizzate la matrice di 
densità di potenze spettrali e la decomposizione ai valori singolari che vengono 


brevemente introdotte in questa sezione. 
La matrice di densità spettrale G,,,(@w) di output assume la seguente forma: 


PSD;1(%) CPSD;z(w) .. CPSD,y(%) 
CPSD;,(%w) PSD,3(%) 


Gyy(®) È: DICCI . aus FICILI 
In cui i termini sulla diagonale principale sono le PSD cioè le funzioni di densità spettrale 
relative a ciascun segnale misurato, mentre i termini fuori dalla diagonale principale sono 


le funzioni di densità spettrale incrociate relative a ciascuna coppia di segnali registrata 


CPSD (Cross Power Spectrum Density). 


Richiamando il Teorema di Wiener che definisce lo spettro di densità di potenza come la 
trasformata di Fourier (FFT Fast Fourier Transform) della sua funzione di 


autocorrelazione si definisce la PSD come: 
PSD(w) = FFT|C,x|? 
In cui la funzione di autocorrelazione C,,(t) è definita nel seguente modo: 
Cax(6) = f x(t)x(t + T)dt 
Mentre la CPSD(w) è definita come: 


CPSD(w) = FFT|Cxy|° 
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In cui la funzione di correlazione incrociata Cxy(1) è definita nel seguente modo: 


Cxy(t) = fox +7)dT 


Per una generica time-history si possono diagrammare le frequenze PSD con i valori 
singolari e si perviene alla realizzazione di un grafico mediante il quale attraverso la lettura 


dei picchi è possibile estrarre le varie forme modali. 


In questo caso studio è stata adottata la decomposizione ai valori singolari detta SVD 


(Single Value Decomposition) della risposta spettrale in cui: 


Gyy(f) = UMISUN' 


In cui la matrice U = [u;, u., ...] è la matrice dei vettori singolari e la matrice [s;] è la 
matrice diagonale dei singoli valori. Come risultato di ciò, andando a graficare i valori 
singolari della matrice di densità spettrale si ottiene un grafico sovrapposto delle densità 
auto spettrali delle coordinate modali. Vale la pena notare che la matrice singolare è una 


funzione delle frequenze. (Brinker, Anderson, & Jacobsen, 2007) 


Attraverso questa procedura, una forma modale viene identificata dove il primo valore 
singolare ha un picco. Nel grafico di Figura 86 si riporta a titolo di esempio un grafico 
dell’SVD. Sono riportati sei valori singolari. Per l’identificazione modale si adotta SV1. Da 
un esame di questa figura si nota che per questo valore sono presenti prevalentemente 


cinque picchi. 
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Figura 86: Grafico SVD di esempio 


3.4.1 Identificazione dinamica 


La tecnica di identificazione dinamica modale Output Only maggiormente impiegata 
oggigiorno è la Stochastic Subspace Identification. Questa metodologia è stata proposta 
nel 1996 da Van Overschee e De Moor che insieme al libro pubblicarono alcuni file Matlab 
per convincere i lettori dell’efficacia dell’applicazione di tale algoritmo come strumento 


forte ed efficiente per l’analisi modale operativa. (Van Overschee & De Moor, 1996) 


Van Overschee e De Moor propongono un modello riportato come quello descritto nella 


Figura 87. 


vi 
System 


ij 


Figura 87: Modello dinamico proposto da Van Overschee & De Moor 
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In cui si indicano con: 
e uxcheèla misura input cui è soggetto il sistema dinamico all’istante k; 
ev, cheè la misura del rumore cui è soggetto il sistema dinamico all'istante K; 
ev, cheèla misura di output del sistema dinamico all'istante k; 

Il sistema strutturale è dotato di massa M, smorzamento C e rigidezza K. 

Sono state individuate due metodologie di risoluzione del problema: 


e Identificazione deterministica del sottospazio in cui sono noti sia l’input uz che 
l'output y,, ed è assente il rumore v,. Attraverso l'applicazione del Teorema 
principale si estraggono gli stati direttamente dai dati di input e di output e quindi 
attraverso l'applicazione di algoritmi si risolvono le matrici incognite del sistema. 

e Identificazione stocastica del sottospazio in cui non è noto l’input (u, = 0), 
quindi si calcolano le matrici del sistema stocastico a partire solamente dai dati in 


uscita. 


3.4.1.1 Identificazione deterministica del sottospazio 


Uno schema concettuale dell’identificazione deterministica del sottospazio è mostrato in 
Figura 88. In questo caso si è considerato il sistema finito lineare e invariante nel tempo. 
| vettori evidenziati da un cerchio sono le misure disponibili mentre v, e wy sono disturbi 
sconosciuti. Il simbolo 4 rappresenta il ritardo. Un'importante osservazione può essere 
fatta in merito all’ipotesi di linearità. Questo approccio descrive bene processi industriali, 
ma in molti altri campi di applicazione, come quello oggetto di questo lavoro di tesi, si 
hanno risposte fortemente non lineari. L'esperienza ha dimostrato che aumentando 
l'ordine del modello di parecchio si riesce a tener conto di fenomeni complessi e si riesce 


ad approssimare bene il problema. 
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Figura 88: Schema di identificazione deterministica del sottospazio 


Il sistema di equazioni che risolve l'algoritmo deterministico è dato dalla seguente coppia 
di equazioni: 
Xk+1 = AxXx + Buy + Wx 
Vk 3 CXx + Dux + Vk 


| vettori: 


e xx è il vettore di stato del processo al generico istante k e contiene i valori 
numerici di n stati; 

e v,ewxsonoisegnali dei vettori non misurabili, si presume che siano sequenze 
vettoriali di rumore bianco a media zero e stazionarie. 


Le matrici: 


e Ab la matrice dinamica del sistema e descrive tutta la dinamica che è stata 


misurata dai sensori. (Brincker & Andersen, 2006) 


0 I 


A3|lmug -M-c 


In cui le componenti sono state introdotte nel precedente paragrafo. 


e Bèlamatrice di input che rappresenta la trasformazione lineare con cui gli input 


influenzano lo stato successivo; 


98 


BS lai 


e Cè la matrice di output che descrive come lo stato interno viene trasferito 


all’esterno nelle misure yy 


e D viene definito come il termine di passaggio. Nei sistemi temporali continui 


questo termine è pari a 0 


Risolvendo il sistema con i dati di input e di output ed applicando i teoremi descritti da 


Van Overschee e De Moor si perviene alle caratteristiche modali del sistema. 


3.4.1.2 Identificazione Stocastica del Sottospazio SSI 


Nell’approccio di identificazione stocastica del sottospazio a differenza del caso 
precedentemente analizzato solo l'output è misurato ed è questo che è stato fatto nel 


presente lavoro di tesi. 


Essendo l’input sconosciuto e se la struttura è soggetta a un rumore bianco la risposta del 


sistema sotto questa sollecitazione conduce alle informazioni modali. 


Uno schema concettuale dell’identificazione stocastica del sottospazio è mostrato in 


Figura 89. 


Figura 89: Schema di identificazione stocastica del sottospazio 


99 


Il sistema di equazioni che risolve l'algoritmo stocastico non tiene più conto del parametro 


di input in quanto per ipotesi u, = 0 e si trasforma quindi in: 
Xx+1 3 Ax + Wx 
Vxk = Cxx + Uk 
I vettori e le matrici sono i medesimi che sono stati descritti al paragrafo precedente. 


Andando a considerare l'equazione dinamica di un sistema generico e sottoposto ad una 


generica forzante F e un generico spostamento y: 
M5(t) + Cy(t) + Ky(t) = F(0) 
Riformulando il modello nello spazio degli stati attraverso la seguente posizione: 


x) 
yo = io 


La prima equazione proposta è definita equazione di stato e modella il comportamento 
dinamico della struttura. La seconda equazione è chiamata equazione di osservazione o 


di uscita. (Mgller, Gade, & Herlufsen, 2005) 
x(t) = Ax(t) + Bu(t) 
y(t) = Cx(t) + Du(t) 


L'ordine della matrice A definisce l’ordine massimo che il modello descrive e questo viene 
scelto in modo arbitrario. La sua scelta determina il massimo numero di modi che il 


sistema può gestire. 


Ai fini dell’implementazione dell'algoritmo sui dati sperimentali acquisiti si riformula il 


modello con riferimento a un dominio tempo-discreto. 
Xx41 = AXx + Wxk 


Vk 3 CXx + Uk 
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Come già accennato, le caratteristiche dinamiche del sistema sono contenute nella 
matrice A di dimensione n x n. Dato un vettore di input w, la matrice trasforma lo stato 
del sistema descritto da n x 1 vettori di stato xx, in un nuovo stato xx +1. La dimensione 
n dello spazio di stato è chiamata dimensione dello spazio di stato. Il vettore di risposta 


del sistema yy è dato dall'unione della parte osservata (matrice C) e dello stato. 


I due vettori che prevedono la risposta del sistema xy e yy, sono correlati attraverso un 


filtro di Kalman per sistemi lineari e invarianti nel tempo: 
Xr+1 = Aîx + Kxek 
ex = Yk — CÎx 
In cui: 
e Kx, èdetta matrice di guadagno di Kalman non stazionario; 


e ex: termine di innovazione modellabile come un processo stocastico del tipo 


rumore bianco gaussiano a media nulla. 


La stima di xy avviene attraverso la tecnica della proiezione ortogonale, maggiori dettagli 


sono stati sviluppati da Van Overschee e De Moor, (Van Overschee & De Moor, 1996). 


Lo step successivo consiste nella risoluzione delle matrici A e C e per le residue sequenze 
Wx e Vx. Infine, per completare il modello equivalente a covarianza completa nel dominio 
tempo-discreto viene stimata la matrice di guadagno di Kalman Ky. Il filtro di Kalman 
presenta un transitorio iniziale la previsione dello stato non è non stazionaria. Quando la 
matrice di stato A è stabile il filtro entrerà in uno stato stazionario e quando questo accade 
la matrice di covarianza del vettore degli stati previsti diventa costante. (Mgller, Gade, & 


Herlufsen, 2005). 
Rielaborando quindi le equazioni si ottiene il sistema dello spazio di stato aggiornato: 
Xr+1 = Aîx + Kxek 


Vk 3 CÈ, + ex 
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Si perviene all'ordine appropriato di un sistema usando la decomposizione del valore 
singolare (SVD Single Value Decomposition) che è una delle più potenti tecniche di 
decomposizione di matrici algebriche lineari. L'ordine del sistema (il numero dei modi) è 
l’unico parametro definito dall'utente nel metodo di identificazione del sottospazio e 
deve essere selezionato con attenzione per ottenere risultati significativi.  (Shokravi, 


Shokravi, Bakhary, Koloor, & Petru, 2020) 


Da un punto di vista teorico ogni valore reale di una matrice mxn può essere 


decomposto come: 
W,0W, = USV® 
In cui: 


e W,eW,sonomatrici di ponderazione che possono essere determinate mediante 


tre metodi: 


* 


«* Principal Component PC; 


* 


«* Unweighted principle component UPC; 
«* Canonical variate analysis CVA; 


e UeV sono matrici unitarie ortogonali; 


e Sèuna matrice diagonale che contiene i valori singolari reali. 


3.4.1.3 Identificazione Stocastica del Sottospazio — Principal Component 


L'algoritmo Principal Component (PC) è uno dei metodi attraverso i quali è possibile 


determinare le matrici di ponderazione. 


Questo algoritmo determina le componenti principali attraverso la matrice di covarianza 
di Toeplitz C; la cui determinazione dettagliata è stata sviluppata da Van Overschee e De 


Moor (Van Overschee & De Moor, 1996) che è data da: 


102 


Da questa formulazione si trova che la decomposizione del valore singolare della matrice 
singolare C;, rivelerà l'ordine del sistema. Nel dettaglio, le colonne hanno dimensione /; 
che è la matrice di estesa osservabilità, mentre le righe 4;“che rappresenta la matrice di 


estesa controllabilità. 


L'algoritmo delle componenti principali considera le matrici di ponderazione W, e W, così 


definite: 
W, = hi 
1 
W,=YT © ? - Y, 
2 P Gul P 


La proiezione delle matrici di ponderazione è espressa come: 


1 


= .& 2 î 
RAI Prrrp] Ta di 


Dove: 
e Yp rappresenta la proiezione ortogonale degli eventi passati 
e Yp rappresenta la proiezione ortogonale degli eventi futuri 


I valori singolari e i vettori singolari di sinistra della proiezione ponderata sono 
esattamente uguali ai valori singolari e ai vettori singolari di sinistra della matrice di 


covarianza di Toeplitz C;. 


3.4.1.4 Diagramma di stabilizzazione 


Nella tecnica SSI uno strumento utile per l'individuazione dei modi propri di una struttura 
è rappresentato dal diagramma di stabilizzazione. Nel dettaglio, attraverso questo grafico 
si può effettuare la separazione tra i modi fisici e i modi spuri nell’algoritmo di 


identificazione. 
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Si definiscono modi stabili (o fisici) quei modi che rimangano costanti al crescere 
dell’ordine del modello. Viceversa, si definiscono modi spuri (o rumorosi) quei valori che 


mostrano tendenze anomale. 


La definizione dell'ordine del modello assume molta importanza. In particolare, una 
crescita dell'ordine del modello conduce all’introduzione di un numero maggiore di modi 
spuri, ma d'altro canto un maggiore ordine del modello permette una più facile 


individuazione degli allineamenti dei modi stabili. 


Si deve scegliere pertanto un adeguato ordine del modello che per il caso studio è stato 


assunto pari a 120. 


La definizione dei modi stabili dipende dalla definizione dei limiti dei seguenti parametri 


modali: 
e Variazione delle frequenze naturali < 1%; 
e Variazione dello smorzamento < 2%; 
e MAC tra due forme modali 99%; 


In riferimento a quest’ultimo punto, si specifica che la MAC (Modal Assurance Criterion) 
è un indicatore statistico ed è calcolato come il prodotto scalare normalizzato tra due 
forme modali consecutive. 

2 
di dil 


MAC= ——————————+6&6 
9-9; ||; - $;l 


Nel dettaglio, questo parametro è stato usato per quantificare l'accuratezza tra le forme 
modali identificate. Valori della MAC prossimi allo zero sono sinonimo di modi ortogonali 
e quindi non c’è corrispondenza tra i due, viceversa quando questo rapporto tende 
all'unità si ha una completa coerenza tra le due forme modali. (Pastor, Binda, & Harcarik, 


2012) 


104 


Una migliore qualità del diagramma di stabilizzazione può condurre all’individuazione di 
allineamenti ben definiti e questo dipende anche dall’algoritmo di identificazione 
utilizzato, dai valori dei parametri di input dell’algoritmo e anche dalla rumorosità delle 


time-history registrate. 
Il diagramma di stabilizzazione è così organizzato: 


e  Nell’asse delle ascisse sono riportati i valori delle frequenze naturali stimate 


dall’algoritmo; 
e  Nell’asse Y, l'ordine del modello. 


Vengono qui riportati i risultati dell’applicazione dell'algoritmo appena descritto su una 
registrazione accelerometrica della durata di 20 minuti ottenuta tramite il sistema di 


monitoraggio della Chiesa di Santa Maria in Via. 


In Figura 90 e in Figura 91 si riportano i diagrammi di stabilizzazione espressi sia in termini 
di frequenze in relazione con l'ordine del modello assunto pari a 50 in questo caso, sia in 
termini di frequenza in relazione con lo smorzamento. Si nota anche che l’algoritmo 
proposto è in grado di cogliere bene frequenze molto vicine tra loro come è possibile 


notare tra 5,5 e 6,5 Hz, in cui si possono distinguere limpidamente tre allineamenti. 
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Diagramma di stabilizzazione 
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Figura 90: Diagramma di stabilizzazione frequenza vs ordine del modello. 
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Figura 91: Diagramma di stabilizzazione frequenza vs smorzamento. 


Si precisa che il diagramma di stabilizzazione di per sé non è sufficiente all’identificazione 
dei parametri modali. Esso rappresenta semplicemente uno strumento grafico per aiutare 
a selezionare le forme modali che più probabilmente rappresentano i modi fisici della 


struttura. (Magalhaes, Cunha, & Caetano, 2009) 
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3.4.2 Estrazione automatica dei parametri modali 


Nel sistema SHM, il ruolo cruciale è svolto dall’automatizzazione del processo di 
identificazione dinamica. In particolare, per ogni registrazione, che come specificato al 
paragrafo 3.3.1 ha durata di dieci minuti, si devono determinare le caratteristiche modali 


secondo il procedimento descritto al paragrafo 3.4.1. 


Ai fini dell’identificazione automatica dei parametri modali della struttura possono essere 
utilizzati algoritmi parametrici che sono focalizzati su tre aspetti complementari: 


(Magalhaes, Cunha, & Caetano, 2009) 


e Concezione di algoritmi di identificazione che possano fornire diagrammi di 


stabilizzazione più chiari; 
e Studio di parametri aggiuntivi per una selezione più accurata dei modi stabili. 


Per tutte queste ragioni, nel presente lavoro di tesi è stata fatta un'analisi cluster per 


interpretare direttamente le informazioni presenti nel diagramma di stabilizzazione. 


3.4.2.1 Il clustering 


L’analisi di cluster (o analisi di gruppi) è una tecnica che si basa sul raggruppamento di 
elementi con caratteristiche simili. Questi gruppi risultanti si definiscono cluster. Questa 
tecnica viene adottata nel contesto dell’interpretazione del diagramma di stabilizzazione 
con l’obiettivo di raggruppare i modi stimati che presentano le medesime caratteristiche 


modali. 


Gli algoritmi di clustering più comuni possono essere divisi in due categorie 


principalmente: 
e Algoritmi di Clustering Non gerarchico; 


e Algoritmi di Clustering Gerarchico. 
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Gli algoritmi di clustering non gerarchico usano solitamente l'algoritmo k-means, in cui il 
set di dati viene diviso in un numero di cluster (K) definito arbitrariamente. Per ogni cluster 
si definisce un centroide cioè un punto al centro del cluster. Quindi l'algoritmo 
inizialmente calcola la distanza euclidea tra ogni dato e ogni centroide. Ogni set di dati 
quindi sarà assegnato al centroide la cui distanza risulti minima. A valle di 
quest’operazione si possono formare nuovi cluster. L'algoritmo procede iterativamente 
finchè i centroidi non si modificano cioè si raggiunge la convergenza per cui non si hanno 


più modifiche nei cluster. 


Il problema maggiore nell’uso di questo tipo di algoritmo consiste nel fatto che si devono 
definire preliminarmente il numero di cluster. Inoltre, questo algoritmo sceglie 


casualmente la distanza dal centro e questo potrebbe condurre a risultati non ripetibili. 


L'algoritmo di clustering gerarchico è basato sulla costruzione di un albero di fusione 


binario che può avvenire attraverso due tecniche: 


e Clustering Gerarchico Agglomerativo (AHC Agglomerative  Hierarchical 


Clustering); 
e Clustering Gerarchico Divisivo (DHC Divisive Hierarchical Clustering) 


Il clustering gerarchico agglomerativo AHC si basa sul fatto che inizialmente ogni dato è 
considerato un cluster. Nello step successivo i due cluster più vicini sono accoppiati in un 
nuovo cluster aggregato. In successione il numero di cluster tende a diminuire fino a che 
tutti gli elementi non sono stati raggruppati e conseguentemente le aggregazioni 


diventano sempre più ampie. (Nielsen, 2016) 


In Figura 92 si riporta un esempio di clustering agglomerativo. 
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Figura 92: Esempio di Clustering agglomerativo — Nielsen. 


Il clustering gerarchico divisivo funziona in maniera inversa. In particolare, parte dal set di 


dati completo e divide progressivamente questo set fino ad arrivare ai singoli elementi. 


La rappresentazione grafica di un’analisi cluster si chiama dendrogramma e in Figura 93 


se ne riporta un esempio, con un'osservazione anche sul tipo di clustering gerarchico. 
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Figura 93: Dendrogramma di esempio. 
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Nel creare un algoritmo di clustering gerarchico agglomerativo è fondamentale definire 


due grandezze: distanza e il linkage. 


II calcolo della distanza, tra i parametri modali del diagramma di stabilizzazione, può 


essere effettuato attraverso le seguenti formule. 
e Distanza euclidea, calcolata come la differenza tra due frequenze f; e f;: 
Deuclidea = [fi = fi| 


e Distanza euclidea normalizzata rispetto alla j-esima frequenza: 


fi Si 
j 


Deuclidea Normalizzata = | f; 


e Distanza di Magalhes, in cui la misura di somiglianza tra i dati viene calcolata 


considerando le frequenze naturali e le forme modali attraverso la MAC. 


fi Si; 


+ (1- MAC) 


Dmagalhes = | f; 

j 
Infatti, se la distanza tra due modi stimati è corta entrambe le stime presentano stesse 
frequenze e stesse forme modali e quindi rappresentano probabilmente lo stesso modo 


fisico della struttura in esame e possono essere inclusi nello stesso cluster. (Magalhaes, 


Cunha, & Caetano, 2009) 


Gli algoritmi di clustering gerarchico si differenziano a seconda di come viene calcolata la 
distanza tra due cluster. Questa particolare distanza viene definita linkage A e può essere 
di diversi tipi: 

e Linkage singolo (SL) in cui si considera la distanza più breve tra due cluster di dati 


Xe Xj 


A(X;,X;) = min D(x;, x%;j) 


xiEXi, xjEXj 
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e Linkage completo (CL) in cui si considera la distanza più lunga tra due cluster di 
dati Xi e Xj 


A(X,,X;) = = max Po x;) 


xiEXi, XjEX 


e Linkage medio tra i gruppi (GAL) in cui si considera la media tra ogni elemento 


dei cluster X; e X; 


A(X;,%;) = dn ) ) DG) 


xieXiXjEXj 


SINGLE LINKAGE COMPLETE LINKAGE GROUP AVERAGE LINKAGE 
SL CL 7 
Suo uo a 
Pa e % CA \ 
/ \ 


Figura 94: Tipologia di Linkage. 


La metodologia proposta per l’identificazione automatica dei parametri modali, si basa 
sul calcolo della distanza attraverso la formula di Magalhes e il linkage tra i cluster di tipo 


singolo. 


Una volta che è stato applicato l'algoritmo di clusterizzazione si effettua una selezione 
delle soluzioni da considerare. Si è definita una soglia minima di valori da prendere in 
esame affinché una frequenza sia considerata stabile. Tale valore è stato assunto pari a 


30. 


Per fare una scelta ancor più accurata dei modi fisici della struttura, sono state poi 
ordinate le frequenze in ordine crescente e sono state considerate quelle la cui differenza 


con il valore precedente sia inferiore a un altro valore soglia. La scelta di questo parametro 
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è fondamentale in quanto tanto più basso è il limite tanto più alto è il numero di cluster 


risultanti e quindi uno stesso modo fisico potrebbe essere diviso in più cluster. È possibile 


comunque usare un valore basso se la qualità del segnale è buona, cioè non presenta 


molta rumorosità. Nel caso in esame, buoni risultati sono stati raggiunti imponendo come 


distanza limite il valore di 0,01. 


Sono state condotte altre operazioni importanti sul segnale acquisito, in dettaglio: 


Sono state rimosse le frequenze con smorzamenti negativi o maggiori di dieci; 


Sono state eliminate le frequenze con valore nullo; 


Non sono state considerate le frequenze con un numero di soluzioni inferiori al 


valore soglia. 


Attraverso l’attuazione di tutte queste operazioni è possibile avere per un generico 


segnale, un diagramma di stabilizzazione depurato dai modi spuri e quindi è possibile 


individuare degli allineamenti che corrispondono ai modi fisici della struttura. 
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Figura 95: Diagramma di stabilizzazione in seguito alla clusterizzazione. 
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Con riferimento al grafico di Figura 95 per questa generica time history sono mostrati i 
cinque gruppi di frequenze “cluster”, individuati a valle della clusterizzazione. In questa 
analisi l'ordine del modello è stato assunto pari a 50. Tuttavia, vale la pena ricordare che 
nei diagrammi di stabilizzazione adottati nella sperimentazione questo parametro è stato 
assunto a 120 per cercare di attuare una selezione più accurata. Le altre probabili forme 


modali non vengono riportate perché hanno percentuali di rilevamento molto basse. 


Tabella 24: Parametri individuati nel diagramma di stabilizzazione di Figura 95. 


Parametri individuati dal diagramma di stabilizzazione 
Numero di Frequenze  |Percentuale di 
poli stabili | della struttura | rilevamento 


2.648 56.00% 


1° Frequenza 


2°Frequenza| 30 | 3.671 60.00% 


3° Frequenza 5.773 40.00% 
4° Frequenza 6.083 54.00% 
5° Frequenza 6.393 54.00% 


MODEL ORDER 50 


28 
30 
20 
27 
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Per ciascun allineamento sono poi selezionate le caratteristiche modali. Questa 
operazione è stata condotta su ognuna delle time history registrate dal sistema di 


acquisizione dei dati. 
3.4.2.2 Eliminazione dei valori anomali 


Per ogni storia temporale vengono selezionati i parametri modali con i criteri definiti 
precedentemente. L'obiettivo di questa seconda fase è quello di selezionare i modi 


individuati preliminarmente depurandoli da eventuali valori anomali. 


Inizialmente, per i parametri modali identificati nella fase precedente, è stata calcolata 
una seconda distanza. Al fine di ridurre i tempi di calcolo è stata fatta una semplificazione 


andando a considerare una distanza euclidea normalizzata: 


fi fi 


Deucl.Norm. = | f; 
J 
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Dopo questa fase preliminare, sulle distanze cosi calcolate, è stata eseguita una nuova 
clusterizzazione al fine di eliminare eventuali tendenze anomale. Per effettuare questa 
operazione è stato scelto un linkage di tipo singolo (SL) tra i vari cluster, e la distanza 


massima tra due frequenze è stata assunta pari a 0,10. 


Una volta applicata tale procedura sono state scelte come frequenze proprie, quelle che 


vengono identificate nel 50% dei casi. 


A valle di questa procedura si ottiene quindi una matrice in cui per ogni modo stabile sono 


riportati i valori delle seguenti caratteristiche modali: 
e Dataeoradiacquisizione; 
e Temperatura interna; 
e Temperatura esterna; 


e Frequenza: 


e Rapporto di Smorzamento £: 


e Valore medio efficace dell’accelerazione agys (Acceleration Root Mean Square) 


che per un segnale campionato è pari a: 


arms 7 


e Accelerazione massima del segnale. 
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4.1 RISULTATI DELLA SPERIMENTAZIONE 


In questo capitolo, si riportano i risultati ottenuti attraverso il software di calcolo Matlab, 
del monitoraggio di Structural Health Monitoring (SHM) effettuato dal 16 novembre al 6 


gennaio 2021 nella Chiesa di Santa Maria in Via. 


Inizialmente vengono riportate le frequenze ottenute a valle dell’identificazione 
automatica dei parametri modali ottenuti dall’algoritmo SSI-PC ed in seguito vengono 
mostrati i risultati derivanti della seconda clusterizzazione. È stata condotta poi una 
statistica descrittiva per i modi di vibrazione ottenuti ed infine, per ogni frequenza sono 


mostrate le eventuali correlazioni con i parametri ambientali. 


4.1. Risultati del monitoraggio 


Nel periodo oggetto di studio, in seguito all'applicazione dell'algoritmo di identificazione 
automatica, sono state processate 3402 time history. Attraverso l'applicazione 
dell'algoritmo SSI-PC come descritto nel paragrafo 3.4.1.4 è stato possibile individuare i 
modi stabili e per ognuno di essi ottenere i parametri modali. Nel grafico riportato in 
Figura 96 si riportano le frequenze ricavate nel periodo dal 16 novembre al 6 gennaio. Si 
può notare un'assenza di dati, segnatamente nel periodo che intercorre tra il 19 e il 21 


dicembre 2020, in quanto il sistema di monitoraggio è stato sospeso per manutenzione. 
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Figura 96: Risultati delle frequenze attraverso l’identificazione automatica. 


È possibile scorgere la presenza costante di alcune frequenze, mentre altre, nonostante 
che siano state processate secondo i criteri del paragrafo 3.4.1.4, non seguono una 


tendenza ben precisa. 


In dettaglio, per poter definire compiutamente il numero di frequenze stabili, è stata 
effettuata una seconda clusterizzazione secondo i criteri riportati al paragrafo 3.4.2.2 al 


fine di eliminare i dati anomali. 


Nella Figura 97, è possibile notare come in seguito alla seconda clusterizzazione il set di 
dati anomali sia notevolmente ridotto e ciò ha reso possibile l'individuazione di otto 
allineamenti corrispondenti ad otto forme modali. Si può notare la presenza costante 
delle prime tre forme modali, mentre le altre non sono rilevate con la stessa continuità. 


Nel paragrafo 4.2 nell'analisi statistica si riportano maggiori dettagli in merito. 
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Figura 97: Frequenze modali ottenute dalla seconda clusterizzazione. 


Per il medesimo periodo di osservazione attraverso le due termocoppie è stato possibile 


registrare le variazioni di temperatura esterna e temperatura interna. 


Per quanto concerne la temperatura esterna, si precisa che questa è stata rilevata come 
anticipato al paragrafo 3.3.1, attraverso una termocoppia posizionata sul presidio di 
sicurezza e quindi si riferisce alla temperatura dell’acciaio dei profili del sistema e non alla 
temperatura dell’aria. In relazione a ciò, si possono notare escursioni termiche massime 


di 15°C, e valori prossimi allo 0 intorno al 24 novembre e tra il 26 e il 28 dicembre. 


La temperatura interna invece ha un calo più o meno costante eccezion fatta peri periodi 
tra il 4 e il 7 dicembre e tra il 22 e il 25 dicembre in cui si nota una leggera crescita in 


concomitanza con all'aumento della temperatura esterna. 
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Figura 98: Temperatura interna rilevata nel periodo di monitoraggio. 
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Figura 99: Temperatura esterna rilevata nel periodo di monitoraggio. 
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4.2. Analisi statistica descrittiva dei modi di vibrazione ottenuti 
È stata condotta un'analisi statistica sui modi di vibrare ottenuti dal sistema di 
monitoraggio. 


Dal precedente capitolo sono stati individuati otto modi di vibrare. Ognuno di essi 


rappresenta un campione di dati da indagare. 
In primo luogo, per ogni set di dati è stata calcolata la frequenza media definita come: 
N 
=1y 
x)== }) Xi 
u N; i 
i=1 
Dove: 
e Nè ilnumero delle osservazioni per ogni campione di dati 
e x; sono le generiche frequenze del campione di dati. 


È stata calcolata poi la varianza per ogni campione. Questo parametro è dato dalla media 


aritmetica dei quadrati degli scarti: 


1 N 
2_ 2 
o°=——_ xi — (x 
v=7) 409) 
i=1 
La varianza da una misura di quanto i dati statistici si discostano dalla media. Più questo 


valore cresce più i dati sono dispersi. 


Quindi è stata calcolata la deviazione standard definita come la radice quadrata della 


varianza: 


N 
1 
o)= pu) Gi - 100)? 
i=1 


Sono stati calcolati anche ulteriori due parametri che sono la Moda e la Mediana. 
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La Moda è definita come il dato che corrisponde alla massima frequenza. Nel caso in cui 


una funzione abbia più di una moda si definisce plurimodale. 


La mediana invece rappresenta il valore assunto dalle unità statistiche che si trovano a 


metà della distribuzione. Per un campione di dati disposti in ordine crescente se: 
e numero di dati è dispari la mediana corrisponde al valore centrale 
e lnumero di dati è pari corrisponde alla media aritmetica dei due dati centrali. 


La distribuzione statistica che meglio si adatta ai valori medi è la legge di Gauss ed è stata 


applicata ai vari set di dati oggetto di studio. 
La densità di probabilità associata alla distribuzione di Gauss è: 


1 Ipo) i 
Pare 


Alla densità di probabilità è anche associata una probabilità di non superamento che è 


definita come l'integrale della densità di probabilità. 


a 1[x-4(x) 
1 GE 


P(x) = [e@e=| E° 0(x) | dx 


In cui u(x) e o(x) rappresentano la media e la deviazione standard precedentemente 


definiti 


Di seguito si riporta per ogni frequenza la curva gaussiana e in dettaglio i valori di media 


e varianza. 


In Figura 100 si riporta la distribuzione normale della prima frequenza e in dettaglio 


vengono riportate la media e la varianza. 
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Figura 100: Distribuzione normale prima frequenza. 


In Figura 101 si riporta la distribuzione normale della seconda frequenza e in dettaglio 


vengono riportate la media e la varianza. 
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Figura 101: Distribuzione normale seconda frequenza. 
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In Figura 102 si riporta la distribuzione normale della terza frequenza e in dettaglio 


vengono riportate la media e la varianza. 
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Figura 102: Distribuzione normale terza frequenza. 


In Figura 103 si riporta la distribuzione normale della quarta frequenza e in dettaglio 


vengono riportate la media e la varianza. 
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Figura 103: Distribuzione normale quarta frequenza. 
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In Figura 104 si riporta la distribuzione normale della quinta frequenza e in dettaglio 


vengono riportate la media e la varianza. 
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Figura 104: Distribuzione normale quinta frequenza. 


In Figura 105 si riporta la distribuzione normale della sesta frequenza e in dettaglio 


vengono riportate la media e la varianza. 
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Figura 105: Distribuzione normale sesta frequenza. 
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In Figura 106, si riportano i parametri statistici della settima frequenza individuata. Si può 
notare come questa frequenza mostri una tendenza anomala. Essendo determinanti i 


valori massimi, una rappresentazione più consona potrebbe essere offerta da una 


distribuzione lognormale. 


Distribuzione normale 7°frequenza 


Numero didati 
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Figura 106: Distribuzione normale settima frequenza. 


In Figura 107 si riporta la distribuzione normale dell’ottava frequenza e in dettaglio 


vengono riportate la media e la varianza. 
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Figura 107: Distribuzione normale ottava frequenza. 


In Tabella 25 sono sintetizzati i parametri statistici delle forme modali individuate. 


Tabella 25: Sintesi dei parametri statistici. 


Numero 
osservazioni 


[n] 


Percentuale 
di rilevazione 


[%] 


Frequenza 
Media 
[Hz] 


Varianza 


[] 


Deviazione 
Standard 
[-] [] [-] 


Moda Mediana 


8° Modo 


2848 83.72 1.64 0.001 0.037 1.588 1.64 
0.024 2.641 2.67 
0.068 3.146 3.27 
0.037 3.689 3.72 
1800 52.91 3.98 0.001 
1784 52.44 4.40 0.002 0.045 4.349 4.40 
1815 53.35 5.82 0.002 


Nell’istogramma di Figura 108, si riportano in forma grafica le percentuali di osservazioni 


delle frequenze individuate. 
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Figura 108: Rappresentazione grafica della percentuale di frequenze rilevate nel periodo oggetto di studio. 


In Tabella 26, si riporta il confronto tra le forme modali identificate nella fase di 


identificazione dinamica globale e quelle ottenute dal monitoraggio. 


Si può notare che il primo modo derivante dall’identificazione dinamica globale non è 
stato rilevato dal monitoraggio, in quanto i sensori sono stati disposti ortogonalmente alla 


facciata, ed essendo questo modo trasversale nell'altra direzione non è stato colto. 


La prima forma modale rilevata invece dal sistema di monitoraggio ha una frequenza di 
1,64 Hz e viene rilevata nell’83,7% dei casi. Non essendo emersa preliminarmente si 
suppone che probabilmente appartenga ad un modo locale della facciata. Per poter 
essere colta il sistema di sensori proposto, dovrà essere implementato con un 


accelerometro posto ortogonalmente a questa parete. 


Le altre cinque forme modali rilevate dall’identificazione dinamica globale invece vengono 
confermate dal sistema di SHM con una differenza tra le frequenze medie minima. Questo 
aspetto è di fondamentale importanza perché conferma l'efficacia del sistema di 


monitoraggio proposto. 
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Infine, sono state rilevate dal sistema SHM altre due forme modali, la settima e l'ottava 
alla frequenza di 5,82 Hz e 6,18 Hz rispettivamente. Non essendo state identificate 


preliminarmente, dovranno essere oggetto di studio per eventuali sviluppi futuri. 


Tabella 26: Confronto tra le forme modali individuate nell’identificazione dinamica globale e nel monitoraggio. 


CONFRONTO FREQUENZE IDENTIFICAZIONE vs MONITORAGGIO 


Frequenze Identificazione 
Forma modale dinamica globale 
[Hz] 


PP W|W[NIN 
vwiolalwualb 
DIANO 
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4,3. Studio dei modi di vibrare con i parametri ambientali 


Una volta che sono state individuate tutte le frequenze sono state messe in relazione con 


alcuni parametri tra cui: 
e Frequenza- Giorni; 
e Frequenza- Ora; 


| parametri sopra indicati sono stati studiati attraverso il diagramma scatola baffi che è 
una rappresentazione adottata in statistica per descrivere la distribuzione di eventi relativi 
ad un campione esaminato, di cui un esempio è riportato in Figura 109. Gli estremi del 
rettangolo in blu corrispondono al primo e al terzo quartile, mentre il valore in rosso 
corrisponde al secondo quartile cioè la mediana. | baffi sono i prolungamenti del 
rettangolo che portano ad individuare due valori il massimo e il minimo. Con una crocetta 
rossa sono riportati invece i valori anomali. Questa tipologia di rappresentazione è molto 


utile per valutare rapidamente i macro-comportamenti di un campione. 


Figura 109: Esempio di diagramma scatola baffi. 
Inoltre, con dei grafici a dispersone sono stati valutati i seguenti andamenti: 
e Frequenza-Temperatura interna; 
e Frequenza- Temperatura Esterna; 
e Frequenza -Smorzamento; 
e Frequenza-Arms; 


e Frequenza- Accelerazione massima. 
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4.3.1. Analisi primo modo 


In relazione alla prima frequenza dal diagramma a) di Figura 110 emerge una significativa 
variabilità delle frequenze rilevate giornalmente con particolare riferimento alla prima 
settimana di osservazione, in particolare tra il 16 e il 23 novembre si assiste ad una 
progressiva crescita della frequenza media per poi oscillare con meno variabilità per le 
successive tre settimane fino al periodo tra il 23 e il 25 dicembre, cui si assiste ad un calo 
evidente, per poi riprendere a salire negli ultimi dieci giorni di osservazione. Viceversa, dal 
diagramma b) di Figura 110 non emerge nessuna variazione sostanziale delle frequenze 


individuate in relazione all’ora. 


Nei parametri confrontati in Figura 111 e in Figura 112 non emerge nessuna sostanziale 


correlazione. 


1° Frequenza - Giorni 


Frequenza [Hz] 
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Figura 110: Diagramma scatole baffi: a) 1° Frequenza — Giorni; b) 1° Frequenza — ora. 
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Figura 111: c) 1° Frequenza - Temperatura Interna; d) 1° Frequenza - Temperatura Esterna. 
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Figura 112: e) 1° Frequenza — Smorzamento; f) 1° Frequenza — Arms; g) 1° Frequenza Accelerazione - 
massima. 
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4.3.2. Analisi secondo modo 


Per quanto concerne il primo grafico di Figura 113 emerge una leggera variabilità delle 
frequenze giornaliere. In dettaglio, si può notare una diminuzione costante della 
frequenza media nella prima settimana di monitoraggio fino al 24 novembre per poi 
ricrescere, mentre una sostanziale deviazione dal comportamento medio si nota anche 
tra il 28 e il 29 dicembre. Dal secondo grafico della medesima figura si nota una piccola 
diminuzione delle frequenze intorno alle ore centrali del giorno per poi aumentare in 


corrispondenza delle ore notturne. 


Nei parametri confrontati in Figura 114 e in Figura 115 non emerge nessuna sostanziale 


correlazione. 
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Figura 113: Diagramma scatole baffi: a) 2° Frequenza — Giorni; b) 2° Frequenza — ora. 
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Figura 114: c) 2° Frequenza - Temperatura Interna; d) 2° Frequenza - Temperatura Esterna. 
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Figura 115: e) 2° Frequenza — Smorzamento; f) 2° Frequenza — Arms; g) 2° Frequenza Accelerazione - 
massima. 
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4.3.3. Analisi terzo modo 


Per quanto concerne il terzo modo, come si può notare dal primo grafico di Figura 116 
emerge una significativa variabilità delle frequenze rilevate giornalmente. Si nota una 
crescita sostanziale della frequenza media nel periodo dal 23 al 25 novembre in antitesi 
con quanto osservato per il secondo modo di vibrare, mentre si nota un decadimento 
della frequenza media dal 23 al 25 dicembre. Invece, dal diagramma b) di Figura 116 non 


emerge nessuna variazione sostanziale delle frequenze con l’ora. 


Nei parametri confrontati in Figura 117 e in Figura 118 non emerge nessuna sostanziale 


correlazione. 
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Figura 116: Diagramma scatole baffi: a) 3° Frequenza — Giorni; b) 3° Frequenza — ora. 


133 


Frequenza [Hz] 


(S) d) 


3°Frequenza - Temperatura intema 3°Frequenza - Temperatura esterna 


Frequenza [Hz] 
W_ wW ww 


sw 


4 6 8 10 12 14 16 
Temperatura interna [C7] Temperatura esterna [C7] 


Figura 117: c) 3° Frequenza - Temperatura Interna; d) 3° Frequenza - Temperatura Esterna. 
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Figura 118: e) 3° Frequenza — Smorzamento; f) 3° Frequenza — Arms; g) 3° Frequenza Accelerazione - 
massima. 
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4.3.4. Analisi quarto modo 


Per quanto riguarda la quarta frequenza, dal primo grafico di Figura 119 si può notare un 
andamento molto simile a quello della seconda frequenza, in particolare si assiste 
inizialmente ad un calo del valore medio della frequenza fino al 24 novembre e sempre in 
analogia con il diagramma di Figura 113 si nota un netto calo della frequenza media tra il 
28 e il 29 dicembre. Si nota una certa similitudine anche per quanto concerne il secondo 
grafico di Figura 119 con una diminuzione delle frequenze nelle ore centrali e un 


conseguente aumento nelle ore notturne. 


Per quanto concerne i diagrammi riportati in Figura 120 e in Figura 121 non sono evidenti 


significative correlazioni. 
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Figura 119: Diagramma scatole baffi: a) 4° Frequenza — Giorni; b) 4° Frequenza — ora. 
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Figura 120: c) 4° Frequenza - Temperatura Interna; d) 4° Frequenza - Temperatura Esterna. 
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Figura 121: e) 4° Frequenza — Smorzamento; f) 4° Frequenza — Arms; g) 4° Frequenza Accelerazione — 
massima. 
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4.3.5. Analisi quinto modo 


Nella quinta frequenza individuata dal sistema di monitoraggio, si possono trarre le 
medesime considerazioni effettuate per il secondo e il quarto modo. Con un decadimento 
nella prima settimana del monitoraggio per poi rimanere più o meno stabile nelle 
successive tre settimane. Anche qui risulta evidente un decadimento nel periodo tra il 28 
e il 29 dicembre per poi tornare ad assumere valori stabili fino al 6 gennaio. Le medesime 
considerazioni possono essere condotte per la il grafico b) in cui appare un lieve calo della 


frequenza nelle ore diurne per poi crescere leggermente nelle ore notturne. 


Per quanto riguarda i diagrammi di Figura 123, Figura 124 come per gli altri modi di 


vibrare, non sono presenti correlazioni significative. 
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Figura 122: Diagramma scatole baffi: a) 5° Frequenza — Giorni; b) 5° Frequenza — ora. 
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Figura 123: c) 5° Frequenza - Temperatura Interna; d) 5° Frequenza - Temperatura Esterna. 
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Figura 124: e) 5° Frequenza — Smorzamento; f) 5° Frequenza — Arms; g) 5° Frequenza Accelerazione — 
massima. 
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4.3.6. Analisi sesto modo 


Anche nel sesto modo, come mostrato nel diagramma di Figura 125, possono essere fatte 
considerazioni analoghe. | valori più alti vengono assunti attorno al 6 dicembre e tra il 24 
e il 25 dicembre. Viceversa, si assiste ad un sostanziale calo durante la prima settimana di 
osservazione e tra il 26 e il 28 dicembre. Da un punto di vista orario, non sono presenti 
fluttuazioni significative, cosi come per quanto riguarda i parametri ambientali, in Figura 


126 e in Figura 127. 
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Figura 125: Diagramma scatole baffi: a) 6° Frequenza — Giorni; b) 6° Frequenza — ora. 
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Figura 127: e) 6° Frequenza — Smorzamento; f) 6° Frequenza — Arms; g) 6° Frequenza Accelerazione — massima. 
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4.3.7. Analisi settimo modo 


Come evidenziato anche nel paragrafo 4.2, la settima frequenza è caratterizzata da una 
tendenza anomala. In particolare, nel grafico a di Figura 128 dal 21 dicembre, si può 
notare un significativo cambio di tendenza rispetto alle prime cinque settimane di 
monitoraggio. Per quanto riguarda il secondo grafico della medesima figura il valore 
medio resta più o meno sempre costante, però è evidente la presenza di un buon numero 
di componenti anomale. AI pari delle altre frequenze rilevate non risultano evidenti 
particolari correlazioni con i parametri ambientali per il periodo oggetto di studio come si 


evince dai diagrammi di Figura 129 e Figura 130. 
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Figura 128: Diagramma scatole baffi: a) 7° Frequenza — Giorni; b) 7° Frequenza — ora. 
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Figura 129: c) 7° Frequenza - Temperatura Interna; d) 7° Frequenza - Temperatura Esterna. 
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Figura 130: e) 7° Frequenza — Smorzamento; f) 7° Frequenza — Arms; g) 7° Frequenza Accelerazione — 
massima. 
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4.3.8. Analisi ottavo modo 


Anche per questa ottava forma modale, vengono confermate le osservazioni fatte in 
precedenza. Si notano i medesimi punti di minimo osservati per la seconda, la quarta la 
quinta e la sesta frequenza. Anche in questo caso nelle ore diurne si assiste ad una leggera 
diminuzione del valore medio della frequenza per poi crescere nelle ore notturne. Per 
quanto riguarda le correlazioni con i parametri ambientali, globalmente non emerge 


nessuna sostanziale dipendenza come si può notare nei grafici di Figura 132 e Figura 133. 
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Figura 131: Diagramma scatole baffi: a) 8° Frequenza — Giorni; b) 8° Frequenza — ora. 
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Figura 132: c) 8° Frequenza - Temperatura Interna; d) 8° Frequenza - Temperatura Esterna. 
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Figura 133: e) 8° Frequenza — Smorzamento; f) 8° Frequenza — Arms; g) 8° Frequenza Accelerazione — 


massima. 


CONCLUSIONI E SVILUPPI FUTURI 


II presente lavoro di tesi ha visto lo studio e l’implementazione di una procedura 
automatica per il monitoraggio in continuo dei parametri modali che caratterizzano il 
comportamento dinamico della chiesa di Santa Maria in Via situata a Camerino. 
Inizialmente, è stata descritta la prima campagna sperimentale intrapresa e che ha visto 
svolgersi prove di vibrazione ambientale con lo scopo di caratterizzare 
approfonditamente il comportamento dinamico dell'intera struttura. Dopo di che, è stato 
affrontato e descritto il sistema di monitoraggio permanente installato nella chiesa. | dati 
registrati da tale sistema hanno fornito la base per lo sviluppo e l’implementazione della 
procedura automatica di estrazione dei parametri modali, parametri su cui attualmente 


si basa il controllo della salute strutturale della chiesa 


A fronte di ciò è stata condotta una campagna di indagine dinamica che ha permesso di 
identificare la dinamica della chiesa nei confronti delle vibrazioni ambientali. In seguito, è 
stato predisposto un sistema di monitoraggio al fine di tenere sotto controllo l'evoluzione 
nel tempo delle caratteristiche dinamiche della struttura, con particolare attenzione alla 


dinamica della facciata. 


Per indentificare i parametri modali, è stata dimostrata l'efficacia dell'algoritmo SSI-PC 
applicato alle time history acquisite che riesce a cogliere frequenze anche molto vicine tra 
loro. Un importante aspetto ha riguardato lo sviluppo di un algoritmo attraverso il 
software di calcolo Matlab per la procedura di identificazione automatica dei parametri 
modali a partire dal diagramma di stabilizzazione. Questo ha permesso di processare in 
poco tempo, e automaticamente, più di 3400 time history acquisite nel periodo dal 16 
novembre al 6 gennaio; non sarebbe stato possibile acquisire un così vasto set di 
parametri se i dati avessero avuto bisogno dell’intervento di un operatore per la selezione 
dei risultati. Da questo punto di vista si è dimostrata l'efficacia dell’applicazione 
dell'analisi statistica di clustering gerarchico che ha consentito di monitorare nel tempo 


otto forme modali. In particolare, è stato possibile monitorare le cinque frequenze 
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individuate durante l’identificazione preliminare ad ottobre e seguire la loro variazione 


nel tempo in termini di frequenze seppur minima. 


Tuttavia, per cogliere il primo modo traslazionale della chiesa è necessario monitorare in 
futuro un ulteriore punto di misura. In generale, questa prima fase di monitoraggio ha 
consentito di gettare le basi per un posizionamento ottimale dei sensori e per la futura 
normalizzazione dei dati ai fini dell’ottenimento di un modello statistico rappresentativo 


delle condizioni attuali della struttura. 


Le frequenze individuate sono state messe in relazione con diversi parametri come lo 
smorzamento, la temperatura interna ed esterna, l'accelerazione massima e 
l'accelerazione media efficace e sembra che il terzo modo torsionale sia il più sensibile ai 
parametri ambientali. Globalmente, per tutti i modi emergono dipendenze con questi 
parametri che è necessario continuare ad indagare in futuro. Sarà quindi opportuno 
monitorare l’edificio anche nel periodo estivo, in quanto per le sette settimane in cui è 
stato in funzione il sistema di monitoraggio la temperatura esterna, registrata con la 
termocoppia sul presidio di sicurezza, si è mantenuta quasi costantemente tra i 4 e i 10 
gradi. Oltre all'introduzione di un numero maggiore di sensori accelerometrici, per cercare 
di capire meglio le dipendenze dalle condizioni ambientali, in futuro si potrà equipaggiare 


il sistema di monitoraggio anche con una stazione meteo. 
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